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ABSTRACT 



A field investigation wae eonducted to determine the impact 
Bf a laFge (0.51+ ha) siibsurface sewage disposal system on groundwater, 

Groundwatar mounding, influenced by the soil conditions, 
hydraulic loading rate, climatic and natural environmental factors, 
can be predicted by either of two mathematical models (Hantush and 
Sykes) . 

Nitrate, chloride, sodium, calcium, sulphate, hardness and 
conductivity were measured to study the magnitude and extent of the 
groundwater contamination, A modified version of the 'Langmuir 
Equation', was used to determine the adsorption capacity of the soil. 
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1.0 INTRODUCTION 

Large leaching beds are being proposed with greater frequency 
by developers as a means of disposal for the effluent from package 
sewage treatment plants serving small subdivisions. 

Large beds offer the economic benefits of scale as compared 
to individual septic tank - leaching bed installations, and the 
operational benefits of year-round use and low energy input compared 
to spray irrigation. 

Possible limitations of effluent discharge into these beds 
may be imposed by excessive groundwater mounding causing surface flood- 
ing (2, 29) and extensive migration of contaminants, especially nitrates, 
into potable water supplies (2, 3, 12, 26). Because of this, there was 
reluctance on the part of regulatory agencies such as the Ministry of 
the Environment and Ministry of Health, to approve such systems. Since 
few field data were available on the performance of large leaching beds, 
there was a need to determine the environmental impact of subsurface 
sewage effluent disposal from these systems, under actual field condi- 
tions. 

The purpose of this study was to determine the environmental 
impact of such large subsurface sewage effluent disposal systems from 
two main points of view. 

1. To determine the effect of hydraulic loading on the level of 
groundwater under the leaching bed. 

2. To determine the impact of subsurface sewage effluent disposal 
on the natural quality of the groundwater. 

In order to achieve these objectives, a suitable site was 
selected adjacent to the Norwood Sewage Treatment Plant in Peterborough 
County, Ontario. A large leaching bed, 8M- m x 64 m in size with 
2,4-40 metres of trenched pipe was installed in the study area in 
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accordance with the local and provincial regulations. The leaching bed 
was dosed with effluent from the Norwood Extended Aeration Sewage Treat- 
ment Plant. 

Fluctuations in the groundwater level and the movement of 
selected contaminants were monitored by 36 piezometers/single level 
sampling wells*'* and 6 multiple level sampling wells installed within 
and outside the perimeter of the leaching bed. Pertinent data were 
collected during a 146 week study period. 

In the analyses of the collected data, two mathematical models 
('Hantush' and 'Sykes') were investigated for predicting the mounding 
of the groundwater under the leaching bed. A modified version of the 
'Langmuir Equation' was used to determine the adsorption capacity of 
the soil. Contaminant concentration plumes were developed to study the 
magnitude and extent of the groundwater contamination. 

After a description of the selected site and effluent disposal 
system, the analyses results of the collected data, based on the study 
objectives, are grouped and presented in two parts: 

Groundwater Mounding; 

Groundwater Contamination. 

The discussion, summary and conclusions on each of the above 
parts are then dealt with separately, followed by the recommendations. 

From this extensive and relatively long-term groundwater 
monitoring program, valuable field data were obtained which provided a 
good insight on the impact of a large subsurface sewage disposal system 
on a groundwater regime . 

''« The terms piezometer and sampling well are used interchangeably 
in this report. 
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2.0 SITE AND DISPOSAL SYSTEM 

2.1 Site Selection 

Seven sites adjacent to existing sewage treatment plants 
were evaluated. The location adjoining the Norwood Sewage Treatment 
Plant was selected on the basis of the availability of an adequate 
area of land for the construction of the disposal system, the quality 
of the plant effluent, preliminary soil conditions, topography, year- 
round accessibility to the site and distance from the Applied Sciences 
Laboratories of the Ministry of the Environment. 

2.2 Site Location 

Norwood is situated on the Trans-Canada Highway, central 
Ontario route, twenty-seven kilometres east of the City of Peterborough 
and ten kilometres north of the Town of Hastings, at the junction of 
Highway 7 and Highway 45. 

The project site was located south of the Sewage Treatment 
Plant in the southwest quadrant of the Town of Norwood and is shown in 
Figure 1 . 

The Norwood Sewage Treatment Plant uses the extended aeration 
process (oxidation ditch) and is equipped with facility for phosphorus 
removal. The design capacity of the plant is 730,000 L/d, with an 
effluent low in Biochemical Oxygen Demand and Suspended Solids and high 
in Nitrates, tjrpical of a small package plant. The treated effluent is 
discharged into the Ouse River, which flows through the Town of Norwood 
on the northeast side of the plant and enters the Trent Canal System at 
the north end of Rice Lake. 
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2.3 Construction of the Disposal System 

The sewage effluent disposal system used in this study, 
comprised mainly : 

(a) A distribution box to give equal flow to all sections of 
the leaching bed. 

(b) Trenched perforated pipe, 2,440 metres in total length, 
for spreading sewage effluent into the soil. 

(c) A storage tank and a double loop syphon system capable 
of feeding a known quantity of sewage effluent to the leaching bed. 

2.3.1 Distribution Box 

The distribution box was constructed in situ with 'ready -mix' 
concrete and was located in the centre of the proposed leaching bed 
at the appropriate elevation. Figure A on plate I shows the inside 
arrangement of the constructed distribution box. The position and 
design of the distribution box is shown in Figures 2 and 3. 

The box was 1.2 m long, 1.07 m wide and 1.2 m deep. The walls 
of the box were 150 mm thick and there was provision for four holes on 
two opposite sides for eight 100 mm diameter headers placed lower thcin 
and at right angle to the axis of 150 iran diameter main feed-pipe. 

A steel baffle plate was secured close to the feed inlet to 
decrease any turbulence and surge during the loading period. Eight 
90° elbows were attached to the header pipes inside the box and cut to 
the same elevation to ensure equal distribution of feed for all the 
headers . Each header provided feed for ten rows of tiles in each of 
the eight bays of the leaching bed. 
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A. DISTRIBUTION BOX CONSTRUCTION 




B. TRENCH CONSTRUCTION 



PLATE I 
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2.3.2 Leaching Bed 

The leaching bed area was 84 m long and 64 m wide and was 
divided into eight equal bays: A, B, C, D, E, F, G, and H. Each bay 
was 19.2 m long, 30.5 m wide and comprised of ten rows 100 mm diameter 
perforated PVC pipes, 30.5 m in length laid parallel 2.1 m apart in 
0.45 m wide trenches, about 0.15 m above the bottom of the trench. 

Each trench bottom was filled with crushed stone or gravel 
around the pipe and backfilled with top soil in accordance with the 
existing Regulation 229 of the Ministry of the Environment. Figure B 
in plate I shows the site during the construction phase. The layout and 
design of the leaching bed is shown in Figure 2, 

The ground elevation ranged from 30.1 m to 29.3 m (Figure 4) 
giving a surface slope of 0.8 ra over the length of the leaching bed. 
Because of this slope, the trenches were 0.51 m deep at the north end 
and 1.36 m deep at the south end of the bed. Each bay of tiles was fed 
by a 100 mm diameter solid pipe from the distribution box. The headers 
receiving feed in all the bays were kept level for an equal distribution 
of effluent. 

The crown elevation of the perforated pipes varied throughout 
the bed from 29.07 m to 28.80 m. The pipes were provided with a 0.10 m 
slope along their 30.5 m length for gravity flow. The end of each row 
was capped and a 13 mm diameter hole was drilled in the centre of the 
cap. A plastic 90*^ elbow was cemented into the hole and a length of 
tygon tubing from the elbow was extended above ground level, as a vent, 
to prevent the development of airlock during the loading cycle. 
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2.3.3 Feed Tank and Syphon System 

The feed tank was an insulated, fiberglass, cylindrical 
shaped tank, 6.1 m in height, 3.2 m in diameter and had a capacity of 
45,000 L. It was positioned inside the Sewage Treatment Plant boundary, 
over a levelled sand pad, close to the final clarifier. A submersible 
pump rated at 205 L/min was connected to a 2M- hour timer, with 15 minute 
increments and used to fill the tank. It took about 210 minutes for the 
submersible pump to fill the feed tank to the point when the syphon was 
activated. Figure A in plate II shows the feed tank at the time of 
installation. Figure 3 shows the layout and design of the feed and 
distribution system. 

In order to feed the large leaching bed, so that 75% of the 
volume of the tiles are filled with effluent (equivalent to about 
15,000 L), a special system was developed which comprised a double 
loop syphon of 150 mm diameter plastic pipe utilizing an initial 5.5 m 
head in the feed tank. The syphon was activated when the liquid in the 
feed tank reached a level of 0.61 m below the top of the tank. 

The pipe carrying the discharge from the tank to the distri- 
bution box was buried in a trench and had a slope of 0.10 m per 30.5 m 
length, for gravity flow. The flow during the first minute was approxi- 
mately 2,400 L/min and decreased to 1,350 L/min after 20 minutes. The 
tank drained 40,900 L in 30 minutes. The submersible pump remained 
activated for 10 minutes after the feed cycle started to ensure complete 
working of the system. 

2.4 Site Exploration 

A detailed exploration of the test site was undertaken 
including: (i) the determination of the soil types and stratigraphy, 
(ii) the flow characteristics of the natural soil deposit. 
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2.4.1 Percolation Test 

The percolation tests ( 'perc tests') were performed at one 
metre depth according to the Ministry of the Environirient standard 
procedure. The test results are presented in Table 1 and the test 
locations are shown in Figure 4. 

TABLE 1 
PERCOLATION TEST RESULTS 



Test 
Location Percolation Time (min/cm) 



1 0.17 

2 0.23 

3 0.25 

4 0.21 

Average 0.22 

The percolation or 't' time varied from 0.17 to 0.25 min/cm 
which indicates a very permeable sandy soil near the ground surface. 
However, results obtained from other soil tests at the site showed 
that the soil became less permeable with increasing depth. It is 
reasonable to expect a larger *t' time in the soil below the 1 m depth, 




OJ 

I 



_ m _ 



2.4.2 Soil Classification and Identification 

A number of boreholes were made at the site to recover soil 
samples for identification and testing. In addition, the soil strati- 
graphy was determined on the basis of the information obtained from the 
boreholes. In Figure 4, the borehole locations and ground surface 
elevation contours are indicated. In the laboratory, grain size distri- 
bution tests were done for most of the soil samples and the Atterberg 
limits (the liquid and plastic limits) were performed on the slightly 
cohesive and cohesive soil samples. 

The classification of the soil samples, the grain size test 
data, the liquid limit and plastic limit results are summarized in 
Table 2. 

2.4.3 Geophysical Exploration and Soil Stratigraphy 
The depth of the bedrock from the ground surface was 

established by resistivity and seismic techniques, which was undertaken 
by the Hydrology and Monitoring Section of the Ministry of the Environ- 
ment. The results are shown in Figure 5. The geophysical test results 
indicated that the bedrock sloped downward approximately in a north- 
westerly direction toward the Ouse River. The depth to bedrock under 
the leaching bed ranged from about 9 m in the southeast corner of the 
bed to 15 m in the northwest comer of the bed and then to about 18 m 
near the river. 

Based on the information obtained during the field investiga- 
tion from the boreholes , the laboratory test data and the results of 
the geophysical exploration, a soil stratigraphy along the centre line 
of the leaching bed, in the north to south direction, is presented in 
Figure 6 . 



TABLE 2 
GRAIN SIZE DISTRIBUTION AND CLASSIFICATION OF SOILS 



Test 










Coefficient of 


Liquid 


Plastic 


Plasticity 




Location 


Elevation 


Sand 


Silt 


Clay 


Uniformity 


Limit 


Limit 


Index 


Soil Type * 






(m) 


% 


% 


% 


(Cu) 


(Wl) 


(Wp) 


dp) 




5 


27 


.1 - 28.8 


41 


54 


5 


9.25 






^ 


Silty Sand 


a 


26 


.7 - 27.0 


7 


84 


9 


8.64 


24.6 


20.5 


4.1 


Silt 


7 


26 


.7 - 27.0 


3 


79 


18 


4.67 


32.5 


18.8 


13.7 


Clayey Silt 


8 




26.8 


5 


86 


9 


8.90 


- 


- 


- 


Silt 






25.3 


25 


75 


- 


6.40 


- 


- 


- 


Sandy Silt 






23.8 


25 


75 


- 


- 


23.6 


19.7 


3.9 


Sandy Silt 






20.7 


- 


85 


15 


- 


24.1 


20.1 


4.0 


Silt 


Northeast 




27.6 


91 


8 





2.31 


- 


- 


« 


Silty Sand 


Corner 




26.7 


- 


- 


- 


- 


- 


- 


~ 


It 11 






23.0 


30 


68 


2 


4.33 


- 


- 


- 


Sandy Silt * 


Northwest 




27.5 


- 


- 


- 


- 


- 


- 


- 


Silty Sand w 


Comer 




26.6 


- 


-* 


- 


- 


- 


- 


— 


M ft . 






25.6 


- 


- 


- 


- 


- 


^ 


- 


It It 






2i+.7 


66 


34 





2.20 


- 


- 


- 


Fine Silty Sand 






23.8 


- 


- 


- 


- 


- 


- 


■m. 


It It If 


Southwest 




28.3 


23 


70 


7 


11.14 


- 


- 


- 


Sandy Silt 


Comer 




















(slightly cohesive) 






27.4 





92 


i 


7.70 


- 


- 


" 


Clayey Silt 






26.4 


- 


- 


- 


- 


- 


- 


* 


Silt (slightly cohesive) 






25.5 





85 


15 


8.85 


- 


- 


- 


Clayey Silt 






24. 6 


(Hole 


terminated, 


, refusal) 


- 


- 


- 


- 


Silt, Sand and Gravel 


Southeast 




28.3 


- 


- 


- 


- 


- 


- 


- 


Silt (slightly cohesive) 


comer 




27.4 


52 


62 


5 


13.3 


■• 


~ 


■" 


Sandy Silt 
(slightly cohesive) 






26.4 


- 


- 


- 


- 


- 


- 


- 


ft II 






25.8 


(Hole 


terminated, 


, refusal) 


- 


- 


- 


- 


Silt with Sand and Gravel 


50 m South 




28.6 


- 


- 


- 


- 


- 


- 


- 


Silt (slightly cohesive) 


of Bed 




27.7 


- 


- 


- 


- 


• 


- 


• 


ff ft M 






26.7 


(Hole 


terminated , 


, refusal) 


- 


- 


- 


- 


Silt (cohesive). Sand 






















and Gravel 



* According to the Unified Soil Classification System. 
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2.M-.M- Permeability Tests 

The permeability (hydraulic conductivity)-'* of the soil was 
an important factor for this study as it affects the groundwater mound- 
ing and the movement of contaminants in the groundwater. Three methods 
were used to measure the permeability of the soils . 

(1) Laboratory 'Constant-Head* Permeability Tests 

(a) For sandy soil samples, the testing technique can be 
found in many soil mechanics text books (23). The collected soil sample 
was disturbed and was compacted in a small container to a density 
similar to the in situ value. 

(b) For the clayey soil samples, the permeability was 
measured in a 'constant-head' permeameter designed by Chan and Kenney 
(11). 'Undisturbed' Shelby tube soil samples were used for the tests. 
The test results are presented in Table 3. It should be pointed out 
that the laboratory permeability values were measured on small soil 
samples, and generally, the laboratory values would be smaller than 
the actual field values for a large volume of soil. 

(2) Field 'Constant-Head' Permeability Tests 

The field permeability apparatus used at the site consisted 
of constant pressure and flow measuring units. The design of the 
apparatus was similar to that described by Wilkinson (37). The constant 
water pressure was maintained by a water reservoir connected to the 
plastic tube of a piezometer which had been driven into the soil. The 
type used for this permeability test was the Geonor Piezometer. Because 
the dimensions of the porous section of the piezometer were small 
(length 0.31 m, diameter 25 mm), the permeability of only a small volume 
of the soil could be measured. 



* The terms 'permeability' and 'hydraulic conductivity' are used 
interchemgeably in this report. 
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TABLE 3 
SOIL PERMEABILITY TEST RESULTS 



Test 
LO'cation 



m. 
11 



Test 
Location 



Iff 



11 



1^ 



m 



Test 
Location 



11 
» 

m 
m 
m 
ii 
ii 

ii. 

m 



ElevatiO'Ti 



28,8 - 27.1 
25.5 
26.4 

25 . 8 



Elevation, 
(m) 



27.5 
26.6 
24.7 

27.6 
26 . 7 

24.9 
28.0' 
27.0 
25.8 
27.3 
26.4 
25.2 
28.3 
27.3 



Elevation 
(m) 



25.6 
25.6 

24 . 9 
24.7 
24,9 
25,3 
26.1 
25.6 
24 . 8 
25,5 
25,4 
25.5 



Lab. 

Permeability ( cm/s ) 



1.6 X 10 

6.0' X 10 

3.0 X 10 

0.3 X 10' 



=4 
-6 
~6 
-6 



Field 
Permeability ( cm/s ) 
'( ' Cone t ant -Head,.'' Te s t ) 



6.6 X 10 

3.6 X 10 
2.8 X 10 
5.1 X 10 
3.4 X 10 
2.1 X 10 
2,0 X 10 

2.0 X 10 
0.8 X 10 

4.1 X 10 
1.0 X 10 

1.7 X 10 
2.3 X 10 
5.0 X 10 



-5 
-5 
-5 

-5 
-5 

-5 
-5 
-5 
-6 
-5 
-5 
-4 
-5 
-5 



Field Permeability (cm/e) 
'(''Falling- Head'' 'Test) '(''Recovery'' Test) 



1.2 X 10 

4.5 X 10 
2.5 X 10 

3.7 X 10 

1.3 X 10 

1.2 X 10 
1.5 X 10 
0.7 X 10 

2.8 X 10 



-4 
-4 
-4 

-4 
-4 
-4 

=4 



0.6 X 10 

12.0 X 10 

7.0 X 10 

8.0 X 10 



3.1 X 10 



-4 
-4 
-4 
=4 
-4 



-4 



0.1 X 10 

0.4 X 10"^ 
0.1 X 10°^ 



1.2 X 10 
1.0 X 10 



=4 



-4 



0.2 X 10 

2.5 X 10 
0.5 X 10 



-4 
=4 



=14, 
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(3) Field Variable-Head Permeability Tests 

The permeability tests were done in piezometers which were 
also used for the regular groundwater level monitoring and sampling. 

Two methods were used to conduct the tests. In the 'falling- 
head' test, water was put in the pipe of a piezometer and the rate of 
the lowering of the water level was recorded. The field permeability 
of the soil can be calculated by the formula suggested by Hvorslev (19). 
In the 'recovery' (or infiltration) test, water was taken out from the 
pipe of a piezometer, and the rate of the rise of the water leve.l in 
the pipe was measured. The equation suggested by Schmid (30) was used 
to calculate the field permeability. Results obtained in a number of 
piezometers are tabulated in Table 3. 

2.4.5 Soil Porosity 

The porosity of the soil was determined at locations 14 to 17 
(Figure 4), by measuring the bulk density of the soil and the soil water 
content. At these locations, the bulk density of the soil was measured 
by the in situ 'sand cone' method in a small excavation, 1 m below the 
ground surface. The porosity of the soil could be computed by assuming 
the specific gravity of the soil particles being equal to 2.67. 

At several other locations, the b-ulk density and the water 
content of the soil samples taken with Shelby tubes were measured and 
the porosity of the soil could then be calculated. 

Table 4 summarizes the test results, which were later used 
to estimate the specific yield of the soil at the site. 
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TABLE 4 
SOIL DENSITY AND POROSITY TEST RESULTS 



Test 
Location 


Elevation 
(m) 


Soil 
Type 


Water 

Content 

% 


Dry 
Density 
(g/cmS) 


% Saturation 


Porosity 


14 


28.33 


Silty Sand 


12.0 


1.52 


42 


0.42 


15 


28.50 


It ti 


6.9 


1.56 


26 


0.42 


16 


29.13 


ir n 


8.4 


1.78 


45 


0.32 


17 


29.00 


ti 1? 


9.8 


1.63 


41 


0.40 


24 


28.06 


Silt 


22.8 


1.48 


76 


0.44 


24 


28.06 


11 


20.9 


1.65 


91 


0.38 


24 


27.56 


M 


22.1 


1.66 


97 


0.37 


24 


27.56 


1? 


21.8 


1.65 


93 


0.38 


24 


26.56 


It 


27.6 


1.52 


100 


0.40 


25 


27.54 


Silty Sand 


11.8 


1.67 


40 


0.44 


25 


27.54 


ti II 


7.5 


1.72 


30 


0.40 


25 


25.74 


Sandy Silt 


18.7 


1.82 


100 


0.34 


25 


25.74 


11 11 


18.2 


1.85 


100 


0.34 
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2.5 Groundwater Monitoring System 

In order to monitor the level of the water table and sub- 
surface movement of the contaminants, a network of piezometers/sampling 
wells was installed inside and around the perimeter of the leaching bed. 
The location and arrangement of the wells were based on the direction 
of the groundwater flow and anticipated movement of the groundwater 
contaminants. Two types of wells were used which are shown in Figure 7. 

Type *A» wells comprised a wellpoint with three rows of 13 mm 
diameter holes down its 0.45 m in length, wrapped with No. 60 mesh 
screen and had a tapered end. Galvanized metal pipe sections, 30 mm 
in diameter and 1.5m long were joined to the well points by galvanized 
couplings to the desired length for the intake of groundwater. 

Type 'B' wells were provided with sample intakes at three 
levels and comprised 0.6 m long, 76 mm diameter PVC pipe sections, 
joined together with plastic couplings. No. 14- rubber stoppers were 
used to separate the sections and provide a watertight seal from each 
other. An end cap sealed the bottom section of the well. A number of 
10 mm diamter holes were drilled into the sides of each section for a 
length of 0.2 m and then covered by No. 270 nylon mesh. A 10 mm 
diameter polyethylene sampling tube was inserted through the side wall 
at the lower end of the section and a 10 mm diameter polyethylene vent 
tube was inserted at the upper end of each sealed section. The top 
section was extended above the ground surface with a sample tube running 
down the centre of the pipe , and being open to the atmosphere did not 
need a vent tube. Figure 'B' in plate II shows an installation of 
Type 'A' (left) and Type 'B' (right) sampling wells. 
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Both types of wells were installed in auger holes which were 
refilled with the native soil. Table 5 shows elevation at the middle 
of the sample intake for the single and multi-level piezometers/ 
sampling wells, in relation to the ground surface and lowest initial 
water table. 

Initially, in June 1976, there were eight wells installed 
within and six wells installed around the perimeter of the leaching bed. 
Later, a network of wells was established on the east side, 10 to 60 
metres away from the bed. By December 1976, three rows at 30, 60 and 90 
metres were added on the west side of the bed. In May 1977, the series 
of wells in the east side had to be removed at the request of the land- 
owner. On May 17, 1977, six multiple level sampling wells (Type 'B') 
were installed adjacent to the existing wells in the bays A, B, D, E, F, 
and G, inside the perimeter of the leaching bed. A total number of 36 
T3^e 'A ' piezometers /wells and 6 Type 'B' piezometers /wells were finally 
used to monitor the level of groundwater table and contaminants movement, 
The network of piezometers and sampling wells is shown in Figure 8. 

Water from the wells was aspirated into one-litre glass 
bottles, using a hand operated vacuum pump. A battery operated depth 
probe was used to monitor the depth of groundwater table in the wells. 

A weekly monitoring and sampling program from July 13 to 
September 28, 1976, established the background data for the project. 
Fluctuations of the water table before subsurface disposal were plotted 
from this background information. The background chemical concentra- 
tion in the groundwater was also established by these data. 
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TABLE 5 
ELEVATIONS OF PIEZOMETER/SAMPLING WELL INTAKES 



Piezometer 


Ground 


Initial 


Location 


Surface 


Water Table 




(m) 


(m) 


1 


29.13 


26.21 


2 


29.46 


26.10 


3 


29.59 


26.00 


4 


29.71 


25.04 


i 


29.15 


26.30 


6 


28.88 


26.27 


7 


29.19 


26.31 


8 


28.17 


25.76 


9 


29.52 


25.83 


10 


29.71 


26.34 


11 


29.66 


26.48 


12 


29.42 


26.48 


13 


29.30 


26.40 


Ik 


28.95 


26.02 


15 


27.77 


25.44 


16 


29.90 


26.00 


17 


30.09 


25.71 


18 


29.70 


26.43 


19 


29.29 


26.35 


20 


29.26 


26.39 


21 


29.21 


26.24 


22 


28.00 


25.84 


23 


30.06 


26.42 


24- 


29.79 


26.55 


25 


29.64 


26.41 


26 


29.54 


26.39 


27 


30.14 


25.55 


28 


30.08 


26.44 


29 


30.14 


26.50 


3Q 


29,98 


26.46 


31 


29.60 


26.37 


32 


29.44 


26.43 


33 


29.50 


26.49 


3k 


29.41 


26.44 


35 


29.11 


26.42 


36 


27.66 


26.33 



Mid-Point Intake Elevation of 
Single Level Piezometer/Well 
(m) 



25.78 
25.75 
26.00 
25.91 
25.74 
25.19 
25.62 
24.41 
25.78 
25.88 
25.97 
25.77 
25.61 
25.23 
24.45 
26.00 
25.64 
26.05 
25.18 
25.60 
25.55 
24.75 
25.36 
25.15 
24.88 
24.79 
25.64 
25.49 
25.45 
25.30 
24.80 
25.11 
25.85 
25.80 
26.06 
25.33 



MID-POINT INTAKE ELEVATIONS OF MULTI-LEVEL PIEZOMETER/SAMPLING WELL 



Ground 
Surface 



Level 
1 



Level 
2 



Level 
3 



It 

i 
c 

B 



29.55 


26.81 


26.18 


25.61 


29.92 


26.82 


26.26 


25.65 


30.07 


27.02 


26.41 


ii* 


30.04 


26.99 


26.36 


25.77 


30.08 


26.93 


26.30 


25.74 


29.85 


28.30 


27.69 


27.11 


29.63 


27.62 


27.04 


26.43 
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2.6 Initial Water Table 

The elevation of the water table, before subsurface disposal, 
was established by nieasuring the depth of the groundwater level in the 
piezometers from the ground surface. The water table contours are 
shown in Figure 8. 

These contours indicate the presence of a slight water table 
mounding below the leaching bed at the beginning of the field testing 
program. On the south side of the bed, a slight depression of the 
water table was noted. This anomaly could have been due to the presence 
of a layer o-F sand, gravel and clayey silt mixture at the piezometer 
tip elevation, and a downward flow of groundwater at this location. 
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3.0 GROUNDWATER MOUNDING 

3.1 Precipitation, Evapotranspiration and 
Water Table Fluctuation 

Groundwater levels fluctuate naturally from season to season 
within the yearly climatic cycle. The geographical location, the 
hydrogeology of the site, climatological influences such as, precipita- 
tion and evapotranspiration all combine to determine the maximum height 
of the water table. Disposal of sewage effluent into the soil through 
a leaching bed is an additional source of influence. 

3.1.1 Precipitation 

Precipitation falls as rain for most of the year and as snow 
or rain mixed with freezing rain during the winter months. The amount 
of precipitation varies from year to year and for the period studied, 
(July 1976 to April 1979), the precipitation total was 2360 mm, giving 
an annual average of 809 mm, at the Norwood site, (13). Although, 
snowfall is measured in centimetres, it is recorded as equivalent milli- 
metres of water. The effect of rainfall on groundwater levels is 
dependent on the intensity and duration of the individual event, but 
the single most significant factor affecting the level of the water 
table is the spring thaw period. In winter, when the mean air tempera- 
ture falls below O^C, moisture and free water freezes in the upper soil 
layers while snow accumulates on the ground surface . As long as the 
frozen layer remains in the soil, the melting surface snow would not 
percolate to the water table. Only at the moment of final melting of 
the frost layer does any free standing surface water seep through the 
soil to the water table. The thawing of the frozen soil may take 
several weeks to complete and by itself, causes an increase in the 
groundwater levels. During a mild early spring, most of the surface 
snow melts quickly while the frost in the soil thaws more slowly. 
This enables much of the snow melt to 'run off before it has a chance 
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to infiltrate into the soil. During a cool late spring, the surface 
snow cover remains longer on the ground allowing the underlying frost 
layer more time to thaw. In other words, the thawing of the frost layer 
occurs at the interface of the frozen upper soil layer and the subsoil 
which is at above freezing temperatures. This produces a gradual upward 
thawing of the frost layer toward the ground surface. Cool climatic 
conditions leave a greater quantity of snow remaining on the ground 
surface and available to percolate to the water table as snow melt when 
the soil is free of frost. Therefore, in general, a cold spring tends 
to give a greater rise of water table than a mild spring. 

In the spring time, the maximum load to the leaching bed was 
equal to the actual hydraulic load plus the quantity of frost and snow 
melt released into the soil from the start of the thaw until the ground- 
water attained its maximum level. The duration of the spring thaw at 
Norwood varied from 4 to 6 weeks. After the spring water level peak, 
the water table declined until the mound reached equilibrium. Minor 
fluctuations in groundwater levels, which occurred throughout the summer 
and fall seasons, were caused by precipitation and evapotranspiration. 

Table 6 shows the pertinent weather data for the winter and 
spring seasons during the study period. As will be seen in subsection 
3.1.3, the maximum rise of the water table during the spring thaw period 
could be correlated to the weather conditions, such as, mean monthly 
temperature, snow cover at the start of the spring thaw, etc. 
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TABLE 6 
PERTINENT WINTER WEATHER DATA AT NORWOOD 





YEAR 


Weather Data 


1976/1977 


1977/1978 


1978/1979 



November (Mean Monthly 

Temperature ( C)) 
December 


- 0.7 

- 9.9 


2.9 

- 7.3 


1.6 
- 4.6 


January 


- 12.8 


- 11.1 


- 9.0 


February 


- 7.3 


- 12.0 


- 13.5 


March 


i.a 


- 5.2 


0.9 


April 


f.l 


3.0 


5.7 


Freezing Degree (°C) day 


908 


1071 


800 


Recorded Winter Snow (cm) 


127.0 


210.0 


145.0 


Snow Cover at Start of Thaw 
(cm) 


i,0 


41.0 


20.0 


3.1.2 Evapotranspiration 







Evapotranspiration (evaporation and transpiration) is a 
difficult parameter to measure accurately and in this study could only 
be estimated. From his study (32), Solomon concluded that the poten- 
tial evapotranspiration is a percentage of pan evaporation. Pan 
evaporation is defined as the gross measurement of evaporation from an 
open pan measured during the frost free period from April to October. 
The pan evaporation values for the Lindsay area Cpart of Simcoe and 
Kawartha lake region) 70 km west of Norwood, for the three year study 
period was 2928 ram of water giving a yearly average of 976 mm (13). 
The mean annual potential evapotranspiration for the Simcoe and 
Kawartha lake region was 580 mm (8), which was approximately 60% of 
the measured pan evaporation. 



TABLE 7 
MEAN MONTHLY TEMPERATURE, PRECIPITATION, AND PAN EVAPORATION AT NORWOOD FROM 1976 TO 1979 (STUDY PERIOD ONLY) 



Month 



January 

February 

March 

April 

May 

June 



Mean Temperature (°C) 
1976 1977 1978 1979 



19.2 



- 9.0 

-13.5 

0.9 

5.7 



July 


19.1 


20.5 


20.2 


August 


18.1 


18.1 


19.1 


September 


13.8 


14.5 


13.4 


October 


5.3 


7.4 


7.4 


November 


- 0.7 


2.9 


1.6 


December 


- 9.9 


- 7.3 


- 4.6 



Precipitation (mm) 
1976 1977 1978 1979 



- 


49.5 


108.7 


- 


40.1 


i.i 


- 


71.6 


56.1 


- 


53.1 


56.0 


- 


57.4 


98.3 


86.9 


58.9 


26.4 


10.5 


70.4 


50.2 


42.4 


124.0 


83.1 


79.8 


102.1 


84.8 


49.0 


65.8 


61.6 


22.1 


114.0 


62.0 


33.5 


114.3 


53.6 



78.5 
58.6 
52.6 
77.0 



Pan Evaporation (mm) 
1976 1977 1978 



63.0 


107.0 


115.0 


73.0 


178.0 


141.0 


164.8 


171.0 


186.0 


173.0 


182.0 


212.0 


150.0 


164.0 


170.0 


165.0 


82.0 


136.0 


64.0 


61.0 


58.0 



1979 



70.5 



03 



- 32 - 

Brandes (7) calculated that the annual evapotranspiration in 
Orillia (same general area as Norwood) was 526 mm. This would indicate 
that the annual evpotranspiration for Norwood is in a range between 
50% - 60% of the annual pan evaporation. If it is assumed that at 
Norwood the evapotranspiration was equal to 55% of the pan evaporation, 
then the amount of evapotranspiration can be estimated. Table 7 
summarizes the monthly temperature, precipitation and pan evaporation 
data for the study period at Norwood. Table 8 shows the daily average 
precipitation and the estimated evapotranspiration. It can be seen that 
at some periods of time the amount of precipitation was exceeded by the 
amount of evapotranspiration. This natural 'deficit' influenced the 
effective net hydraulic loading rate in the leaching bed which in turn 
caused minor fluctuations in the water table . 



TABLE 8 
COMPARISON OF PRECIPITATION WITH ESTIMATED EVAPOTRANSPIRATION 



YEAR 




















Precipi- 
tation 

mm 


Est'd 
Evapotran- 
spiration 
/day 


Precipi- Est'd 
tation Evapotran- 
spiration 
mm /day 


Precipi- Est'd 
tation Evapotran- 
spiration 
mm/day 


April 


1.04 


2.10 


(d)A 


1.77 


1.96 


(d) 


1.96 


2.11 (d) 


May 


3.30 


2.37 


(s) 


1.85 


3.61 


(d) 


3.17 


2.50 (s) 


June 


2.90 


3.02 


(d) 


1.96 


3.13 


(d) 


0.88 


3.41 (d) 


July 


3.56 


3.07 


(s) 


2.27 


3.23 


(d) 


1.62 


3.76 (d) 


Aug. 


1.37 


2.66 


(d) 


4.00 


2.91 


(s) 


2.68 


3.02 (d) 


Sept. 


2.66 


3.03 


(d) 


3.40 


1.50 


(s) 


2.83 


2.49 (s) 


Oct. 


1.58 


1.14 


(s) 


2.12 


1.12 


(s) 


1.98 


1.03 (s) 




* (d) = 
(s) = 


deficit 
surplus 
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3.1.3 Water Table Fluctuation 

Typical water table fluctuations are presented in Figure 9. 
Two groups of results are presented: one from piezometers installed 
inside the leaching bed and the other from piezometers installed at 
varying distances from the bed. The water table below the leaching bed 
was influenced by the hydraulic loading from the leaching bed plus the 
natural environmental factors such as, precipitation and evapotranspira- 
tion. The water table fluctuations at a certain distance from the centre 
of the bed were less influenced by the hydraulic loading. To explain 
the groiindwater fluctuation data, the results for piezometers 30 and 2 
are shown in Figure 10. Piezometer 30 (solid line) was located within 
the leaching bed area while piezometer 2 (dotted line) was situated about 
90 m from the perimeter of the bed and least influenced by the hydraulic 
loading. In Figure 10, the significant water table fluctuations are 
numbered for piezometer 30 and the following notes with the correspond- 
ing numbers outline some details and explanation of the fluctuations. 

1. Week to Week 12 

A weekly monitoring program was initiated on July 14, 1976. 
The natural water table declined slowly after the construction of the 
leaching bed. The depth of the water table below the ground surface at 
week 11 (elevation = 26.46 m) was used as the base value for the calcu- 
lation of the subsequent net groundwater levels. 

2. Week 12 to Week 17 

Hydraulic loading (Feed) to the leaching bed was started on 

October 5, 1976, at a rate of 122,700 L/d which is equivalent to 

2 . 
22.8 L/m of l(?aching bed area. Within the first week of loading, there 

was a significant increase in the groundwater levels beneath the bed, 

although, this increase was not even throughout the tiled area. At the 
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south end of the bed, the water table rise was 1.64 m in one week then 
continued increasing slowly over the next four weeks by an additional 
0.46 m. The water table rose at a more even rate at the north end of 
the leaching bed, increasing by 0.43 m, 0.49 m, 0.46 m, 0.49 m, and 
0.2 m during this five week loading cycle. 

3. November 10, 1976 

On this date , liquid was observed bubbling through the top 
soil and ponding on the ground surface above the bed. The system was 
hydraulically overloaded (i.e. in a state of failure) and loading was 
terminated. (See plate III). The average water table rise below the 
centre of the bed was 1.99 m. The total volume of feed injected into 
the leaching bed was 4.296 x 10° litres. 

4. Week 17 to Week 34 

The artificial groundwater mound declined between November 
10, 1976 and March 8, 1977. This was a rest period (i.e. no feed). 
The response of the water table at the south end of the bed was initial- 
ly different from that at the north end. The water level fell . 82 m 
in the south and 0.55 m in the north during the first week of the rest 
period then continued to decline slowly until the water table was at 
pre-injection levels. 

5. Week 34 to Week 40 

A sharp rise in the water table was noted between March 8, 
1977 and April 19, 1977, although, there was no feed to the leaching 
bed. This was the natural increase in water levels caused by the spring 
thaw period and was general over all the monitored area. The average 
increase caused by spring thaw was 0.76 m. 
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A. BED FAILURE 




B. BED FAILURE 



PLATE III 
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6» Week UP to Week Ul 

Feed to the leachi.ng bed was re-started on April 21, 1977, 
at 40,900 L/d, Unfortunately, an electrical malfunction developed 
in the timer allowing the feed pump to operate continuously. A total 
of 1.06 X 10 litres of secondary effluent was injected in a four-day 
period. The response of the water table to the massive load was an 
increase of 1.37 m at the south end of the bed and 0.64 m rise at the 
north end. The average increase was 0.98 m within the bed and 0.52 m 
at the perimeter. There was very little or no increase away from the 
tiled area because of the short duration of the loading period. 

7. Week 41 to Week 44 

This was a short rest period to allow the groundwater mound 
to decline. 

8. Week 44 to Week 53 

Hydraulic loading was resumed on May 17, 1977, at a rate of 
40,900 L/d. The groundwater level was 1.01 m above the base eleva- 
tion at the bed centre and a subsequent rise in the water table was 
noted below the leaching bed area. 

9. Week 53 to Week 54 

The system was shutdown for one week on July 19, 1977, at 
the request of the Water Resources Branch of the Ministry of the 
Environment which was conducting a river survey downstream from the 
treatment plant. Consequently, the water table declined at this time. 
10. Week 54 to Week 73 

The feed was continued at 40,900 L/d from July 25, 1977, 
to December 9, 1977, with an increase in the water table. 



- 39 - 



11. Week 73 to Week 74 

An ice blockage in the line from the pump to the feed tank 
caused the system to be shutdown for one week on December 9, 1977, 
while repairs and modifications were made. The groundwater mound 
declined. 

12. Week 74 to Week 87 

Loading was resumed December 16, 1977. The feed was main- 
tained at 40,900 L/d with a corresponding mound in the groundwater 
levels. From week 44 to week 87, including 2 one-week shutdown periods, 
the water table rose 1.86 m below the bed centre. The overall trend 
from week 44 through 87 was an increase in the height of the water table. 
A feed rate of 40,900 L/d was adequate to sustain an average 1.42 m rise 
above the base groundwater level, (i.e. elevation = 26.46 m). Volume of 
feed applied to the system from week 44 to 87 was 12.025 x 10 litres. 

13. Week 87 to Week 91 

On March 14, 1978, the water levels increased over all the 
monitored area due to the spring thaw. The artificial mound under the 
leaching bed disappeared as the spring thaw reached its peak. Exclud- 
ing the centre of the bed which already contained the artificial mound, 
the average increase in the water table during the spring thaw was 1.69 m. 
Between April 12 and April 18, 1976, water was observed bubbling up 
through the soil onto the ground surface. (See plate III). However, 
the feed was continued since it was felt that the flooded condition 
would last for only a very short time. When the groundwater levels 
peaked on April 18, 1978, the average height above the base level over 
all the monitored area was 2.73 m. The rise of the water table during 
the spring thaw in 1978 was greater than that in 1977. This observa- 
tion is in agreement with the weather data in Table 6. In 1978, the 
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freezing degree days, the recorded winter snow and the amount of snow 
cover at the start of the thaw were all greater than those in 1977. 
14. Week 91 to Week 100 

This was a period of decline and stabilization of the water 
table after the spring peak levels had occurred. The feed to the 
leaching bed was maintained at 40,900 L/d. The decrease over all 
the field excluding the centre of the leaching bed during this 10 week 
period was 1.4U m. Within the area of the bed the decline in the water 
table was only 0.98 m, due mainly to the continuing influence of the 
hydraulic loading. On May 16, 1978, liquid was again noticed breaking 
out from the ground surface. (There had been a prolonged period of 
steady rainfall at this time combined with subsoil still partly saturat- 
ed from the spring thaw). The surface area remained very moist for 
some time after the flooding and breakout of liquid. The rise in ground- 
water levels from the base value on June 14, 1978, was 1.77 m at the bed 
centre. Total volume applied to system from week 87 to week 100 was 
3.722 X 10^ litres. 

15. Week 100 to Week 126 

This was a continuous loading period of 40,900 L/d from 
June 14, 1978, until the final shutdown of the system on December 13, 
1978. The groundwater levels below the leaching bed remained relatively 
steady with minor fluctuations due to sporadic rainfall and evapotran- 
spiration. The water table rise above the base elevation (i.e. 26.46 m) 
at week 126 was 1.54 metres at the centre of the leaching bed. Total 
volume of liquid injected into the bed for 26 weeks was 7.157 x 10^ 
litres. 
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16. Week 126 to Week 138 

The feed to the leaching bed was stopped on December 13, 1978. 
The groiindwater levels immediately began to decrease and continued to 
drop throughout the winter until the spring thaw of 1979. 

17. Week 138 to Week 144 

This was the annual spring thaw period. The natural water 
levels increased during six weeks and peaked on April 18, 1979, at 
1.54 m above the base elevation. The average increase caused by the 
thaw period over all the monitored area was 1.15 m. The water table 
rise for the 1979 spring thaw was smaller than that in 1978 as the 
recorded winter snow and the snow cover was less in 1979. (See Table 
6). The natural water levels were declining again by week 146 of the 
study. 

In Table 9, a summary of the net fluctuations of the ground- 
water table with respect to the initial base elevation is presented. The 
quantity of liquid applied to the leaching bed throughout the study 
totalled 28.26 x 10^ litres. 

3.2 Mathematical Techniques for Predicting 
Mounding 

A number of mathematical models in various forms are avail- 
able today for the prediction of the rise and fall of groundwater mounds 
under recharge areas subject to percolation including some early 
research (5, 15, 17). In the last decade, additional models were 
developed (18, 24, 25, 29, 31, 36). 

Generally, in the models, an aquifer system is represented by 
a set of partial differential equations. With given governing boundary 
conditions and aquifer parameters, the differential equations can be 
solved both exactly and numerically. In the Norwood Project, Hantush's 
closed-formed solutions and Sykes ' finite-element numerical model were 



TABLE 9 

SUMMARY OF GROUNDWATER FLUCTUATIONS (IN METRES) IN DIFFERENT GROUPS OF 
TEST LOCATIONS AT VARYING DISTANCES FROM THE LEACHING BED 



Week No. 



- 


- 12 


12 - 


- 17 


17 - 


- 34 


34 - 


- 40 


40 - 


- 41 


41 - 


- 44 


44 - 


- 87 


87 - 


- 91 



91 - 100 

100 - 126 

126 - 138 
138 - 144 



Comment 



Pre start-up (Base) 

High loading rate 

Rest period 

Spring thaw 1977 

Overload 1 week 

Rest period 

Low loading rate 

Spring thaw 1978 
(+ low load rate) 

Decline (stabilize) 
(+ low load rate) 

Summer period 
(low loading rate) 

Rest period 

Spring thaw 1979 



GROUP 1 

(center) 

Net Cumula- 
tive 



0.00 0.00 

+1.99 1.99 

-2.04 -0.05 

+0.72 0.67 

+0.98 1.65 

-0.64 1.01 

+0.85 1.86 

+0.89 2.75 

-0.98 1.77 

-0.23 1.54 

-1.16 0.38 

+1.11 1.49 



GROUP 2 
(Perimeter) 



Net 



Cumula- 
tive 






00 





.00 


+1 


48 


1 


48 


-1. 


44 





14 


+0 


73 





77 


+0 


52 


1 


29 


-4 


6 





83 


+0 


50 


1 


.33 


+1 


67 


3 


00 



-1.58 1.42 



-0.42 1.00 



-0.56 
+1.14 



0.44 
1.58 



GROUP 3 
(30 metres) 



Net 



Cumula- 
tive 






00 





.00 


+0 


84 





84 


-0 


80 





04 


+ 


70 





74 


+0, 


09 





83 


-0 


20 





63 


+ 


48 


1 


11 


+1. 


63 


2 


74 



-1.44 



1.30 



-0.41 0.89 

-0.51 0.38 
+1.12 1.50 



GROUP 4 
(60 metres) 



Net 



Cumula- 
tive 



0.00 





00 


+ 0.40 





40 


-0.40 





00 


+0.81 





81 


+0.12 





83 


-0.19 





64 


+0.23 





87 


+1.83 


2 


70 



-1.50 1.20 

-0.57 0.63 

-0.26 0.37 

+1.17 1.54 



GROUP 5 
(90 metres) 



Net 



Cumula- 
tive 



0.00 0.00 



+0.86 0.86 

-0.01 0.85 

-0.18 0.67 

+0.15 0.82 

+1.62 2.44 

-1.23 1.21 

-0.65 0.56 

-0.20 0.36 

+1.22 1.58 



m 
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used to analyse the field data. The main reason for selecting the 
Hantush model was the relative simplicity of its solutions which would 
make it easier to use for more routine design problems. Sykes' model 
is flexible in many ways. This model represents the more recent 
numerical approach to the groundwater mounding problem. 

In this section, some details of Hantush 's solutions are 
presented and in the appendix, the mathematical treatment of Sykes' 
model is given. 

The equations developed by Hantush can be used to estimate 
the rise and fall of a groundwater mound in response to a uniform 
vertically downward rate of percolation for both rectangular and 
circular recharging areas. The solutions were developed from the 
Laplace equation in terras of the head averaged over the depth of 
saturation. 

The following simplifying assumptions are made in the model :- 

1. The aquifer is homogeneous, isotropic, infinite in areal 
extent and resting on a horizontal impermeable base. 

2. The formation coefficients and the rate of percolation are 
constant in time and space. 

3. The constant rate of deep percolation, relative to the 
horizontal flux, is so small that the groundwater flow lines are 
almost completely refracted in the direction of the tilt of the water 
table. 

4-. The water table is initially horizontal and remains below 
the bottom of the recharge area. 

5. The relative rise of the water table to the initial water 
table depth is less than 50%. 

6. The Dupuit-Forschheimer assumptions are valid. 
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The flow problem can be approximated by the following 
boundary-value problem, (see Figure 11 for co-ordinate system). 



— ^ + -^ + i-r;) f(x, y) 
3X^ 3Y^ ^ 



W^ at 



subject to the following initial and boundary conditions: 
Z(x, y, 0) =0 



8Z(0, y, t) 
9x 



3Z(°°, y, t) 

3x 



W = 



W < X < a, 
X > a, 



az(x, 0, t) = 
By 

9Z(x, °°, t) = 

ay 

0<y<b, t>0 
y > b, t > ... 



where 



V = 

b = 

Z = 

k = 

s = 

W = 

t = 

a = 

b = 



kb/s 

0.5 |h, (0) + h(t, 

constant of linearization 



.5 fh^ (0) + h(t^)] 



(1) 



(la) 
(lb) 

(Ic) 

(Id) 
(le) 



permeability (hydraulic conductivity) 

specific yield 

constant rate of percolation 

time since incidence of percolation at water table 

dimension of basin in x direction 

dimension of basin in y direction 
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FIGURE 11 

Geometry and Symbols for Rectangular and Circular Infiltration 
Areas and Underlying Groundwater Mound in Unconfined Aquifer (17) 
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Successive applications of Laplace and Fourier transforms 
and the respective inversion formulas yields the following for the 
first quartile. Due to symmetry, this result is applicable to other 



quartiles . 



Z(x, y, t) = ^ 



erf (^^ ^ erf (^^] ^ ^ 

\/577/ \Av t/J 



where A = 2w/k 



(2) 



After changing the variable of integration from x to ti the 
rise of the water table can be expressed as: 



2 2 

h - h. = wvt 

^ 2k 



f ^' (^^ ^) t s. (^^, ^ 



(3) 



where D = /4vt 



s''- (<=^, 3) = 



erf (- 



a 



/ 



-) . erf (■ 



^ 



-) dx 



Values of the above function are tabulated by Hantush. The 
corresponding equation for the decline of the water table if percola- 
tion ceases is: 



h^ - h^2 , 2(x, y, t) - Z(x, y, t - t^) 



(4) 



where t^ is the time when water ceased to arrive at the water table 
and Z(x, y, t) and Z(x, y, t - t^) are the right-hand side of equation 
(3) with t and t - t as time factors, respectively. This solution of 
hypothetically superposing on the flow system a discharge equal to the 
uniform rate of percolation at t = t^ is possible since the differen- 
tial equation is linear in h^. 
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For a circular recharge basin, Hantush based his solution on 
a Laplace equation of the form, 

3^Z ^ i az ^ 2W .f(r) = 19Z (5) 

9j^2 r 3r k r 9t 

with the following initial and boundary conditions: 

Z(r, 0) =0 (5a) 

9Z(0, t) = (5b) 



= (5c) 



W = 





dr 












9Z( 


°°, 


t) 












9r 






W 







< 


r 


< 


R 











r 


> 


R 



(5d) 
(5e) 



where R = radius of circular recharge area. 

The problem was solved using the Laplace and zero order 
Hankel transformations to produce the final equation. 

h^ - h^^ = (2V/7Tk) f (q, p) (6) 

2 
where V = wit R 

q _ vt 

r2 

p _ r^ 
R 

As for the case of a rectangular recharge area, the decay of 

the water table for a circular area after cessation of infiltration is 

similarly given by, 

h^ - h^^ = Z(r, t) - Z(r, t - t^) (7) 

Hantush 's equation for the rise of the mound below a 
rectangular recharge basin can be written in a simpler form (6) and 
is easier to use than Hantush 's equation in its original form. 
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The equation for the rise of the mound in uncon fined 
aquifers below rectangular basins is: 

F { (W/2 + x)n, (L/2 + y)n } 

+ F { (W/2 + x)n, (L/2 - y)n } 

+ F { (W/2 - x)n, (L/2 + y)n } 

+ F { (W/2 - x)n, (L/2 - y)n }{ (8) 



-H \^ / 



where h = height of water table above impermeable layer at 
X, y, t 

X, y and time t (see Figure 12) 

H = original height of water table above impermeable layer 

V = arrival rate at water table of water from infiltration 
a 

basin 

t = time since start of recharge 

f = finable porosity (1 > f > 0) = specific yield = s 

L = length of recharge basin (in y direction) 

W = width of recharge basin (in x direction) 

n = (i+.t.T/f)"'^ 
1 
F(a, 3) = J erf (at ) . erf (Bt )dT {which was tabulated by Hantush} 
& 

The equation for the decay of the groundwater mound after 
cessation of infiltration is, 

h^ ^^ .-H = Z(x, y, t) - Z(x, y, t - t ) (9) 

X, y, t s 

where t = time since water began to arrive at the water table 

t = time since water ceased to arrive at the water table 
s 

Z(x, y, t) and Z(x, y, t - t ) = right-hand part of equation (8) 
with t and (t - t ) as time factors. 
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3.3 Analysis of Groundwater Mounding 
by Hantush's Method 

In this section, the selection of the values of the soil 
(or aquifer) parameters used in the calculation of the magnitude of 
the groundwater mounding by Hantush's equations will be briefly 
discussed. The calculated results will then be compared with the 
field experimental data obtained in the periods when the hydraulic 
loading rates were 122,700 L/d and 4-0,900 L/d. 

3.3.1 Selection of the Soil Parameter Values 
As noted in 3.2, the values of the fillable porosity, 'f, 
(specific yield, 's'), the permeability, 'k', (hydraulic conductivity) 
and the thickness of the aquifer, 'H', are required in Hantush's 
equations to compute the amount of the groundwater mounding beneath a 
recharge area. These values were measured in the field exploration 
program which was described in detail in Chapter 2, Because the natural 
soil conditions were quite variable and the thickness of the aquifer 
was not constant, some engineering judgement was required in the selec- 
tion of the values of the parameters used in the calculations, 

(1) Fillable Porosity, 'f, (Specific Yield, 's') 

A value of 0.2 was chosen for the fillable porosity, 'f. The 
selection was based on literature data which are presented in Figure 13. 
In addition, the results of measurements on the porosity and the degree 
of saturation of the soil samples (see Table 4) obtained at the test 
site were used as a guide. 

(2) Permeability, 'k^, (Hydraulic Conductivity) 

The test results are summarized in Table 3. It can be noted 
that the test results vary significantly, depending on the testing 
methods and the soil types. As shown in Figure 6, the sandy silt on 
the north side of the leaching bed is more permeable than the clayey 



!>a 




0.001 
Clay 



0.01 



Silt 



1 0.1 1 

I Sand 

Median gram £i2e 



10 
Gravel 



100 
Cobbles 



FIGURE 13 



RELATION BETWEEN MEDIAN GRAIN SIZE AND 
WATER-STORAGE PROPERTIES OF ALLUVIUM FROM • 
LARGE VALLEYS. (REF. P. 394 IN 'HYDROLOGY' 
BY STANLEY N. DAVIS AND ROGER J. M. DeWIEST, 
JOHN WILEY & SONS, INC. 1966) 
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silt on the south side of the leaching bed. Table 10' is a summary of 
the averaged ''k^ re:s.ults for the sandy silt material obtained by 
various types of field tests. 



TABLE 10 
Averaged ''k' Results for Sandy Silt 
k (cm/s) Type of Field Test 

2.3 X 10'"'^ 'Falling head'' 

-4 
5,1 X 10 'Recovery'' 

-4. 
3.7 X 10 'Falling head' and 'Recovery' 

~5 
3 '5 X 10 ''Constant head' 



It should be pointed out that in a near-horinzontal ground- 
water flow from, a gTQimdwater mound at the test site, t'he more permeable 
soil (i.e. the sandy silt material), under the leaching bed would have 
more important influence on the mo^unding. Therefore, the 'k' results 
of the more permeable soil were used in the calculations of the mO'Unding. 
( 3 ) Thickness of the Aquifer 'H' 

It is noted in Figure 5 that the thickness of the aquifer, 'H', 
increases in the direction towards the river. 'However, only a constant 
'H' can be used in Hantush's equations. Therefore, the value of 'H' was 
assumed to be equal to 12 m. This assumption was made from geop'hyslcal 
data recorded in the area from the leaching bed to t'he Ouse River. The 
thickness of the aquifer was equal to the difference in elevations of 
t'he bedrock and the water table at the outset of t'he study. In the 
groundwater mounding calculations, a value of 8 m was also used and this 
value was a]>proximately equal to the thickness O'f the sandy silt beneath 
and in the vicinity of the leaching bed'...,^. 
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3.3.2 Comparison of Experimental and Theoretical 
Groundwater Mounding (Hydraulic Load = 
122,700 L/day) and Decline 

With the values of 'k', 'H' and 'f* (or 's') presented and 
discussed in subsection 3.3.1, the magnitude of the groundwater mound- 
ing was computed for several combinations of permeability ('k') and 
aquifer thickness ('H'). The calculated results were compared with 
the experimental data of three groups of piezometers. The first group 
of piezometers are: 24, 25, 29 and 30, which were located within the 
leaching bed and were closest to the distribution box. The second 
group of piezometers are : 26 and 31 , which were located on the north 
side of the distribution box, and the third group are: 23 and 28, 
which were located on the south side of the distribution box (see 
Figure 8). In the calculations, the field-determined 'k' values (i.e. 
3.7 X 10~ cm/s and 3.5 x 10~ cm/s) were used, and in addition, a 
'back-calculated' 'k' value was also used. This 'back-calculated' 
permeability value was obtained by comparing the field data at piezo- 
meters 24, 25, 29 and 30 with the theoretical results for a number of 

'k' values. By the trial -and-error procedure, a 'back -calculated' 'k' 

-4 . . 

value of about 5 x 10 cm/s was obtained for the case m which 

H = 12 m and f = 0.2. This 'back -calculated' 'k' value was later used 
with the field-determined 'k' values to predict the long-term ground- 
water mound under the bed. 

The calculated results for a number of combinations are 
summarized in Tables 11a and b and are presented in Figures 14a, b and 
c for comparisons with the field data. As shown by Figure 14a, the 

average data of piezometers 24, 25, 29 and 30 compare quite well with 

-4 
the theoretical data for the 'k' value of 3.7 x 10 cm/s. However, 

_5 
the mounding for the 'k' value of 3.5 x 10 cm/s is greater than the 

actual measurement . 



■r 
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TABLE 11a 



COMPARISON OF HANTUSH CALCULATED GROUNDWATER 
MOUND WITH FIELD DATA 



Loading Rate = 122,700 L/d: f = 0.2: H=12Tn 

Location From Bed Centre ~ 16.0 m East or West; 10.7 m North or South 





FIELD DATA (m) 
Average of 24, 
25, 29, 30 






'k^ 




Days 


(m/d) 
( cm/s ) 


5.0 X 10 


0-32 _4 
3.7 X 10 


0.03 
3.5 X 10 


3 




0.38 




0.34 


0.34 


0.35 


7 




0.84 




0.73 


0.73 


0.80 


10 




1.10 




0.98 


1.00 


1.15 


14 




1.28 




1.27 


1.31 


1.60 


16 




1.46 




1.40 


1.46 


1.83 


21 




1.70 




1.68 


1.78 


2.39 


24 




1.78 




1.83 


1.95 


2.72 



Loading Rate = 122,700 L/d: f=0.2: H=8m 

Location From Bed Centre = 16.0 m East or West; 10.7 m North or South 



FIELD DATA (m) 
Days Average of 24, 
25, 29, 30 



(m/d) 
(cm/s) 



0.43 
5.0 X 10 



-4 



0.32 
3.7 X 10 



-4 



0.03 
3.5 X 10 



-5 



3 


0.38 


7 


0.84 


10 


1.10 


14 


1.28 


16 


1.46 


21 


1.70 


24 


1.78 



0.34 


0.34 


0.35 


0.77 


0.76 


0.81 


1.05 


1.05 


1.15 


1.38 


1.40 


1.61 


1.53 


1.57 


1.84 


1.87 


1.95 


2.40 


2.06 


2.10 


2.74 
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TABLE lib 

COMPARISON OF HANTUSH CALCULATED GROUNDWATER 
MOUND WITH FIELD DATA 



Loading Rate = 122,700 L/d: f = 0.2: H = 12 m 

Location From Bed Centre = 16.0 m East or West; 31.0 m North or South 



FIELD DATA (m) 
Days Average of 
26, 31 



(m/d) 
(cm/s ) 



0.H3 
5.0 X 10 



-4 



0.32 
3.7 X 10 



-*+ 



0,03 
3.5 X 10 



-5 



FIELD DATA (m) 
Average of 
23, 28 



3 


0.13 


7 


0.U9 


10 


0.73 


14 


0.96 


16 


1.18 


21 


1.51 


24 


1.55 



0.20 


0.20 


0.25 


1.23 


0.43 


0.44 


0.54 


1.54 


0.59 


0.61 


0.74 


1.64 


0.78 


0.82 


1.00 


1.65 


0.88 


0.95 


1.13 


1.70 


1.09 


1.15 


1.45 


1.76 


1.21 


1.29 


1.64 


1.79 



Loading Rate = 122,700 L/d: f = 0.2: H = 8 m 

Location From Bed Centre = 16.0 m East or West; 31.0 m North or South 



FIELD DATE (m) 
Days Average of 

26, 31 



(m/d) 
(cm/s ) 



0.43 
5.0 X 10 



-4 



0.32 
3.7 X 10 



-4 



0.03 
3.5 X 10 



-5 



FIELD DATA (m) 
Average of 
23, 28 



3 


0.13 


■? 


0.49 


10 


0.73 


14 


0.96 


li 


1.18 


21 


1.51 


m 


1.65 



0.20 


0.21 


0.27 


1.23 


0.44 


0.45 


0.56 


1.54 


0.61 


0.63 


0.77 


1.64 


0.83 


0.85 


1.04 


1.65 


0.93 


0.96 


1.17 


1.70 


1.17 


1.22 


1.50 


1.76 


1.30 


1.37 


1.69 


1.79 
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In Figure 14b ^ the comp'arison of the experimental curve with 
the theoretical curveiS is reasonably good even through the 'k'' values 
are quite different. The better comparison could be attributed to the 
short period of the hydraulic loading and the distance of the piezo'- 
meters from the centre of the bed, (i.e. smaller mound at larger 
distance frorai the centre of the bed). For' piezometers 23 and 28,, the 
theoretical and experiinental results compare poorly (see Figure 14c). 
This co'uld be due tO',: 

(i) different sod,l conditions on the sO'Uth part and the north 
part of the leaching bed, and 

(il) the 'k' values were for the soil below the north part of 
the leaching bed. 

In Tables 11a and b, the results for H = 12 m and 8 m are 
bo'fh presented. As would be expected, groundwater mound would be 
higher if the ''H' is 8 m instead of 12 m,, ''k' and 'f values being 
equal . 

Figures 15a,, b and c, sbow co'iiparisons of the experimental 
and theoretical data for the decline of the water table after the 

hydraulic loading of 122 ,,700 L/d was stopped. In the calculations , 'k'' 

-5 -4 -4 

values of 3.5 x 10' cm/s , 3-7 x 10 cm,/s and 5 x 10 cm/s were used 

to obtain the theoretical results, 'H'' and 'f being 12 m and 0,2, 
respectively. For different permeability values,, the initial elevation 
of the water table before the decline would be different because a 

smaller 'k' value would yield a higher groundwater moiund. Therefore,,, 
in comparing the theoretical and experim,ental data in Figures 15a,, b and 
c, only the rate of the decline is of interest. TwO' co'iiments may be 

made on the results in Figure 15: 

-4 - 5 

(i) For the ranges of the 'k' values (10 to 10 ' cm/s) used, the 

rate in the lowering of the ground,water mound doe,s not seem to differ 
significantly, and 



' .'^ 
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(ii) The measured rate of the decline was somewhat faster than 
what is predicted. 

3.3.3 Groundwater Mounding During Week ^H to Week 126 
During this period, with the exception of a few short durations, 
the hydraulic loading rate was 40,900 L/d . The fluctuations of the 
water table beneath and adjacent to the leaching bed are presented in 
Figure 9. As explained in section 3.1, the movement of the water table 
beneath the bed was influenced by two factors: 

(i) the hydraulic loading from the bed, and 

(ii) precipitation and evapotranspiration at the site. Because 
Hantush's method does not take into account the effect of evapotranspira- 
tion and it also assumes a constant hydraulic loading rate during the 
loading period, it becomes quite difficult to compare the short-term 
minor fluctuations of the water table with the theoretical curve. 
Fortunately, in design, it is the maximum rise of the water table which 
is of practical significance. Therefore, the comparison of the calcu- 
lated maximum rise of the groundwater mound with what was measured in 
the field would be of prime interest. The theoretical rise of the water 
table at piezometer 30 was calculated by Hantush's equation for several 
cases and the results are tabulated in Table 12. According to the theory, 
the initial rise of the water table is fast at the beginning of the 
hydraulic loading then the rise slows down with time. Referring to the 
data of piezometer 30 in Figure 10, the water table elevation 26.46 m at 
week 12 could be considered as the 'base line' for the measurement of the 
net groundwater mound. Therefore, in April 1978, the maximum increase of 
the water table was 2.75 m. This large increase was the result of the 
combined effect of the hydraulic loading to the leaching bed and the 
infiltration of surface water from the 1978 spring thaw. From week 100 
to week 126, the net rise of the water table was maintained at 1.5 m and 
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TABLE 12 

COMPARISON OF HANTUSH CALCULATED GROUNDWATER 
MOUND WITH FIELD DATA 



Loading Rate = 40,900 L/d: f = 0.2: H = 12 m 

Location From Bed Centre = 16.0 m East or West; 10.7 m North or South 





FIELD DATA (m) 


Days 


Average of 23, 




24, 25, 26, 28, 




29, 30, 31 


100 


For 26 weeks 


200 


(week 100 to 


300 


week 126 of 


400 


study) 


500 


1.54 



(m/d) 
(cm/s) 



0.43 
5.0 X 10 



-4 



0.32 
3.7 X 10 



-4 



0.03 
3.5 X 10 



-5 



1.30 


1.53 


3.43 


1.67 


2.00 


5.77 


1.89 


2.29 


7.44 


2.04 


2.49 


8.71 


2.17 


2.64 


9.71 



Loading Rate = 40,900 L/d: f=0.2: H=8Tn 

Location From Bed Centre = 16.0 m East or West; 10.7 m North or South 





FIELD DATA (m) 


Days 


Average of 23 




24, 25, 26, 28, 




29, 30, 31 


100 


For 26 weeks 


200 


(week 100 to 


300 


week 126 of 


400 


Study) 


500 


1.54 



(m/d) 
( cm/s ) 



0.43 
5.0 X 10 



-4 



0.32 
3.7 X 10 



_4 



0.03 
3.5 X 10 



-5 



1.59 


1.84 


3.59 


2.09 


2.47 


6.23 


2.39 


2.85 


8.16 


2.59 


3.12 


9.62 


2.79 


3.32 


10.78 
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There was very little fluctuation in the elevation of the water table. 
However J for piezometer 2, the elevation of the water table during the 
same period was approximately at the same elevation before the hydraulic 

loading of 40,90i0 L/d was applied to the bed. As pointed out before, 
piezometer 2 was not affected significantly by the hydraulic loading in 
the bed. That is to say, the fluctuation of the water table at this 
location was mainly due to precipitation and evapotranspiration. Because 
piezometers 30i and 2 were located not more than 150 m apart, it may be 
ass'umed that the climatic influence on the water table at these locations 
was about equal. Therefore, if the water table at piezometer 30' had not 
been affected by the many weeks of hydraulic loading, the water table 
would have returned close to elevation 2'6.46 m. However, the net rise 
of the water level in piezometer 30' with respect to elevation 2'6.'46 m was, 
1.5M- ra. Tiifrrefores it can be inferred that this groundwater mound was 
the result -of many weeks of hydraulic loading in the leaching bed. 

'The measuped net rise (1.54 m) in the water in day SO'Oi was 
smaller than the predicted values s.ummarized in Table 12. The main 
reasons could be : 

(i) the average effective hydraulic loading was less than 
40,900 L/day; 

'(ii) the field-determined and the ''back-calculated' 'k* values 
were too' small , and 

(iii) the effective thickness of the aquifer was .greater than the 
ass'cuned 12 m. 
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3.4 Groundwater Mounding Analysis by 
Sykes' Finite -Element Model 

The extent of the groundwater mounding under the leaching bed 
at the site was also calculated by Sykes' coinputerized model, which is 
described in the appendix. In using the finite -element computer model, 
the test site area is first divided into many small areas (elements) 
each of which has four comers (nodes) (Figure 16). To facilitate the 
comparison of the computer and experimental results, some of the nodes 
were placed at the piezometer locations. 

In Figures 17 and 18, the computed results by Sykes ' model for 
two loading rates are compared with those of Hantush for a constant 
H = 12 m and the field data. 

Because Sykes' finite -element model has the flexibility of 
varying the input data (e.g. 'f, *k', 'H' and hydraulic loading rate, 
etc.) in each element, it was decided not to use a constant aquifer 
thickness in the additional analyses. Instead, the thickness of the 
aquifer was calculated by subtracting the elevation of the water table 
prior to the testing program from the elevation of the bedrock as 
determined by geophysical techniques. In Figures 19 and 20, the computed 
results are compared with results obtained by the Hantush analyses, in 
which the thickness of the aquifer was assumed to be 12 m. 

As can be seen in Figures 17, 18, 19 and 20, good comparison 
of results was observed between Sykes' and the Hantush model. 

3. 5 Discussion on Groundwater Mounding 

From the extensive groundwater monitoring study conducted at 
Norwood, Ontario, experience and knowledge can be gained in two areas: 
(i) the fluctuations of the w.Ttor tahlo boneath a ] arp.e leachinj\ 
bed including the relative important influence of precipitation and 
evaporation on groundwater mounding, and 
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Figure 16 ■ Finite Elements and Nodes Used In ttie Sykes Model for the Study Areo 



- 63 - 



4.0 , 



FIGURE 17 



< K-5. 0«10'''*on./S f»B. 20, H~12in > 




X DATA AVG. 24, 2S. 20. 30 
HANTUSH 



20 

TIME (DAYS) 

SYKES' COMPUTED RESULTS COMPARED TO HANTUSH 
(LOADING RATE 122700 L/day) 





^> o . 


( k-5. 0»10'"'*oin/S f-0. 20. 


H-12i« 


> 




















3.0 . 










\y 












LU 

q: 












UJ 

-I 
m 


2.0 . 










1.0 . 












< 


._ 


^ „• ■• 




UJ 

< 


















HANTUSH 






0.0 . 








SYKES 






, 


PHMn 


, 


-f 







FIGURE 18 



100 200 

TIME (DAYS) 



300 



SYKES' COMPUTED RESULTS COMPARED TO HANTUSH 
(LOADING RATE 40900 L/day) 



- 64 - 



4. , 



v5 3. .. 

UJ 
O) 

oc 

UJ 



CD 

< 



2. .. 



UJ 

S 1.0 



0.0 







< k-5. B«10~'*o«i/S f-0.20. H-VARIABLE THICKNESS > 






X 


DATA AVG. 

HANTUSH 

SYKES 


24. 


25. 


29. 


30 



10 



20 
TIME CDAYS> 



30 



40 



FIGURE 19 SYKES' COMPUTED RESULTS COMPARED TO HANTUSH 

CLOADING RATE 122700 L/doy> 



4. , 







< 


K-5. e'le^^om/S f-0.20p H-VARIABLE THICKNESS ) 


E 


3.0 . 






UJ 




* 




UJ 

ffi 


2.0 , 
1.0 . 


■ 








UJ 

1- 
< 








HAWTt lew 






0.0 . 


1/ 




SYKES 













100 200 

TIME (DAYS) 



300 



FIGURE 20 SYKES' COMPUTED RESULTS COMPARED TO HANTUSH 

(LOADING RATE 40900 L/day> 



- 65 - 

(ii) the practical use of two mathematical tools, (i.e. Hantush 
and Sykes ) for the prediction of groundwater mounding beneath a 
recharge area or basin. 

As discussed in sections 3.3 and 3.4, both the Hantush and 
Sykes methods gave reasonable, although conservative, predictions of 
the groundwater mounding under a large leaching bed provided that, 
(i) the 'back-calculated' 'k' value or the right selection of 
the field-determined 'k' value was used, and 

(ii) for the Hantush model, the rapid groundwater rise caused by 
a large quantity of infiltration at the time of spring thaw was not 
considered because of the model's requirement of a constant hydraulic 
loading rate . 

These assumptions lead to a number of practical considerations 
in the prediction of the groundwater mounding under a recharge basin 
(or area) such as a leaching bed. 

(i) In the calculation of the magnitude of the groundwater mound, 
values of soil (or aquifer) parameters are required. In this study, 
these values were obtained either by 'back-calculation' from the field 
data or by actual measurements in a number of boreholes made at the site 
The mounding prediction was reasonable when the 'back-calculated' 'k' 
value or the 'k' value determined by the 'falling head' and the 
'recovery' methods was used. As a result of this study, it is felt that 
the 'falling head' and the 'recovery' 'k' tests in piezometers are 
suitable for sandy silt soils. However, it is not certain whether these 
field tests would yield reasonable 'k' results for clayey soils for the 
groundwater mounding calculations. Perhaps, a model-scale field test or 
a field pumping test would give more meaningful 'k' and specific yield 
(s) results for fine-grained, cohesive soil deposits. 
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(ii) In this study, it was found necessary to consider the water 
table rise caused by the spring infiltration in order not to under- 
estimate the maximuni mounding. As indicated in Figure 10, the rise 
of the water table at various seasons was different. This was only 
reasonable because of the stochastic nature of climatic conditions. 
Therefore, in the design, it would be necessary to study weather data 
(precipitation, evaporation, temperature, etc.) at a potential site 
for as many years as practically feasible. In addition, pr^ior to the 
installation of a leaching bed at a potential site, it is useful to 
install a number of piezometers at the site to monitor the natural 
fluctuations of the water table for a minimum of one year. The one- 
year water table fluctuation data, together with the weather data from 
past years can be used to make a reasonable estimate of the maximum 
rise of the water table at the potential site. A statistical approach 
may be used to predict the probability of the occurrence of the maximum 
rise in the water table due to spring thaw infiltration and the calcu- 
lated mound caused by the design hydraulic loading in the leaching bed. 
In selecting the probable maximum rise of the water table due to spring 
infiltration, the designer should consider the chance for the occurrence 
of this maximum and the effect of surface flooding due to the excessive 
rise of the water table . 

As a precautionary measure, a few piezometers may be installed 
near the centre of the leaching bed and the water table during the 
operation of the bed will be monitored. It may also be feasible to 
adjust the hydraulic loading rate in the leaching bed at the time of 
the spring thaw in order to avoid flooding of the bed. 

(iii) In this study, both the Hantush and Sykes methods were used 
for the calculation of the groundwater mound. Because of the simplicity 
of the Hantush method, it would be a more practical tool for routine 
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calculations. Furthennore , at the present time, the inaccuracy in 
the prediction of mounding lies more in the inability to practically 
assess the overall soil properties. 

However, it should be pointed out that the Sykes model has 
a number of flexible features which Hantush's method does not have. 
For example, in each element of Sykes' model, the soil properties, 
the infiltration rate and the thickness of the aquifer can be different, 
Therefore, Sykes' model can be used as a refined design tool if the 
soil properties at a potential site can be measured with reasonable 
accuracy. In addition, the Sykes finite -element model can more 
accurately represent the groundwater loading caused by seasonally 
varying infiltrations. For the Norwood site, the effect of spring 
infiltration was greater than that of the hydraulic loading of 
40,900 L/d. As stated in the discussion of the Hantush model, the 
model assumes zero infiltration outside of the leaching bed. This 
assumption is not necessary with the finite -element analyses. There- 
fore, for sites in which natural infiltration or soil heterogenetics 
are important, a finite -element type analysis may be warranted. 



GROUNDWATER CONTAMINATION 
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4.0 GROUNDWATER CONTAMINATION 

The magnitude and extent of groundwater contamination due 
to subsurface sewage effluent disposal were investigated by conducting 
a regular and comprehensive sampling program to monitor the variation 
of contaminants concentration with time, distance and depth, under and 
in the vicinity of the leaching bed. 

Type 'A' wells were used to monitor the contaminants concen- 
tration inside and outside the boundary of the leaching bed, while the 
multi -level Type 'B' wells were used to measure the concentrations at 
different depths within the perimeter of the leaching bed. 

Since chemical pollutants travel farther and persist longer 
than bacterial pollutants (10), it was decided to restrict this investi- 
gation to chemical contamination of the groundwater. Because of the 
ease of monitoring and typical behaviour characteristics of certain 
chemicals to enter and move freely within an aquifer (2, 20, 21, 27, 
28), this study is mainly focused on subsurface movement of chloride, 
nitrate, sulphate, sodium, calcium, hardness and conductivity. Concen- 
tration levels of other contaminants and parameters in the case of 
sewage effluent feed are included. 

Chemical analyses of all the groundwater samples collected 
were performed using standard analytical methods (16, 33), in the 
Central Laboratories of the Ministry of the Environment, in Toronto. 

All the statistical values and graphs that are presented were 
analyzed on an electronic computer (HP - 982 5A) and are based on 
standard mathematical equations and statistical curves. 
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4.1 Variation in Chemical Concentration 
of Groundwater 

In order to study the effect of the subsurface sewage effluent 
disposal on the quality of groundwater, variation in the chemical concen- 
trations of the groundwater were monitored before and after the disposal 
of sewage effluent. 

4.1.1 Background Concentration (C^) 

A 12 week sampling program of the groundwater was conducted in 
the project area to establish the natural quality of the groundwater 
before subsurface disposal. Because of a wide difference in the back- 
ground concentration values between the north and south ends of the 
leaching bed, it was decided to use the lowest of the background values 
recorded at the individual sampling well locations for calcium, 
sulphate, hardness and conductivity. The field averages were used for 
chloride, sodium and nitrate. The analyses results of these samples are 
presented in Table 13. 

4.1.2 Feed Con centration (Cp) 

The chemical concentration of the feed (before subsurface 
disposal) during the total time of the study, could be grouped under 
three main periods. 

Period 1, covered the time between October 5, 1976 and 
November 10, 1976 when the hydraulic loading rate was 122,700 L/d. 

Period 2, was a rest period between November 10, 1976 and 
May 17, 1977 when the operation of the system was stopped due to bed 
failure. However, the quality of the plant effluent was monitored 
during this time. 

Period 3 , began on May 17, 1977 when the operation of the 
system was re -started and the loading was maintained at a reduced rate 
of 40,900 L/d to December 13, 1978. 
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TABLE 13 



BACKGROUND CHEMICAL CONCENTRATION IN THE GROUNDWATER 
BEFORE SUBSURFACE DISPOSAL 



Average Background Concentrations in mg/L 



Well 
No. 



Nitrate 



Chloride 



Sodium 



Calcium 



Sulphate 



Avg. 0.23 4.1 5.5 
»'f Conductivity reported in ymhos/cm 



36 



16 



Conduc- 
tivity * 



215 



Hardness 



1 


0.10 


2.1 


1.2 


14 


5 


89 


38 


2 


0.11 


2.7 


1.9 


11 


7 


144 


74 


3 


0.17 


4.7 


»+.o 


40 


19 


200 


111 


i+ 


0.70 


4.4 


3.8 


28 


17 


230 


101 


5 


0.10' 


4.4 


2.8 


14 


5 


57 


23 


6 


0.32 


4.2 


1.6 


11 


4 


61 


25 


7 


0.10 


3.3 


2.5 


14 


11 


82 


3S 


8 


0.15 


4.4 


2.0 


15 


11 


85 


35 


9 


0.40 


4.3 


2,7 


29 


9 


160 


88 


10 


0.18 


4.5 


6.2 


38 


14 


208 


114 


11 


0.10 


5.8 


5.7 


44 


25 


280 


147 


12 


0.14 


4.4 


i*.8 


44 


13 


256 


160 


13 


0.10 


4.3 


tf.2 


25 


17 


220 


100 


14 


0.14 


5.0 


1.5 


19 


9 


75 


55 


15 


0.40 


4.6 


if. 4 


17 


11 


117 


46 


16 


0.73 


5.8 


U.2 


62 


23 


335 


194 


17 


0.18 


3.7 


6.5 


103 


20 


578 


366 


18 


0.10 


4.6 


4.8 


49 


20 


307 


165 


19 


0.10 


6.2 


2.8 


42 


19 


255 


135 


20 


0.13 


7.3 


2.3 


25 


20 


220 


120 


21 


0.10 


7.1 


1.9 


25 


17 


215 


100 


22 


0.40 


5.8 


6.5 


20 


16 


133 


54 


23 


0.13 


2.8 


6.1 


120 


30 


564 


328 


24 


0.10 


3.1 


4.7 


32 


16 


209 


107 


25 


0.43 


3.9 


2.9 


35 


21 


208 


109 


26 


0.18 


3.8 


2.8 


45 


25 


265 


144 


27 


0.10 


5.5 


6.8 


89 


18 


438 


258 


28 


0.10 


3.0 


7.8 


94 


21 


549 


306 


29 


0.10 


3.3 


10.0 


91 


26 


527 


297 


30 


0.15 


5.9 


7.1 


67 


26 


392 


210 


31 


0.12 


5.6 


4.8 


39 


21 


240 


124 


32 


0.99 


4.0 


3.7 


23 


15 


164 


78 


33 


0.10 


2.8 


4.5 


27 


12 


195 


103 


34 


0.17 


5.5 


3.9 


28 


16 


195 


90 


35 


0.10 


4.8 


3.3 


32 


11 


355 


185 


36 


0.30 


3.0 


2.7 


26 


7 


210 


81 



110 



TABLE 14 
CHEMICAL CONCENTRATIONS IN THE FEED BEFORE SUBSURFACE DISPOSAL 







PERIOD 1 






PERIOD 3 












Cheniicals and 
Other Parameters " 




Loading Rate 
122 700 L/d 






Loading Rate 
40 900 L/d 






Total for 

Loading Periods 

1 and 3- 








No. of 








No. of 








No. of 










Samples 




Range 


Mean 


Samples 




Range 


Mean 


Samples 


Range 


Mean 




B0D5 


6 




4-12 


6.5 


60 




1-30 


8.7 


If 


1-30 


8.5 




Suspended Solids 


i 




5-15 


6.8 


58 




<5-20 


8.5 


m 


<5-20 


8.3 




Total Nitrogen (as N) 


i 




24-36 


27.4 


56 




13-32 


20.1 


it 


13-36 


20.8 




Nitrate (as N) 


& 




22-32 


25.0 


56 




4-26 


15.5 


m 


4-32 


16.4 




Sodiuin (as Na) 


i 




67-81 


74 


55 




60-100 


75 


ii 


60-100 


75 




Potassium (as K) 


i 




11-15 


12.2 


5if 




6-15 


9.7 


10 


6-15 


10.0 


1 


Chloride (as CI) 


i 




81-104 


90 


53 




75-154 


104 


5i 


75-154 


103 




Calcium (as Ca) 


i 




100-116 


105 


53 




101-139 


116 


59 ■ 


100-139 


115 




Sulphate (as SO4) 


s 




59-79 


68 


56 




41-98 


64 


62 


41-98 


65 




Conductivity (pmhos/cm) 


i 




850-1010 


947 


51 




910-1120 


1009 


57 


850-1120 


1002 




Hardness (as 03003) 


i 




295-306 


297 


55 




280-379 


320 


61 


280-379 


318 




Alkalinity (as CaCoa) 


s 




163-205 


190 


59 




162-300 


226 


65 


162-300 


223 




Magnesium (as Mg) 


2 




9-12 


10.5 


53 




5-11 


7.4 


55 


5-12 


7.5 




p« 


i 




7.3-7.9 


7.6 


58 




6.4-8.6 


7.7 


64 


6.4-8.6 


7.7 





•" All concentrations are reported in mg/L except Conductivity and pH values 
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Table 14 shows the range and mean values of the chemical 
concentrations and other parameters in the feed for periods 1 and 
3 and for the total period, (period 1 and 3), of subsurface effluent 
disposal. 

The variation in the quality of the Norwood Sev/age Treatment 
Plant secondary effluent for the total study period is shown in Table 
15 which indicates the concentration levels of the selected contaminants 
and parameters for 15%, 50% and 85% of the time. 

TABLE 15 

VARIATION IN QUALITY OF PLANT SECONDARY EFFLUENT 
FOR THE TOTAL STUDY PERIOD 



Chemicals and 
Other Parameters 
(mg/L) 



No. of 
Samples 



Values % of Time Equal to or Less Than 



IS% 



8S% 



Total Nitrogen 


m 


16.9 


20.8 


25.7 


Sodium 




77 


66.8 


74.4 


82.9 


Chloride 




m 


81.2 


99.7 


122.3 


Calcium 




m 


102.6 


112.5 


123.3 


Sulphate 




8t 


54.4 


63.8 


74.9 


Conductivi 


ty '^ 


fi 


945.7 


974.6 


1004.3 


Hardness 




83 


282.0 


311.5 


344.2 



* Conductivity reported in ymhos/cm 
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4.1.3 Chemical Concentration of Groundwater After 
Disposal (C) 

Table 16 shows the maximum concentration level of the selected 
contaminants recorded at the sampling well locations in the project site. 
The highest of the maximum recorded concentrations of nitrate, chloride 
and sodium occurred at well location 25; of calcium, sulphate and hard- 
ness at well location 11 and conductivity at well location 27. The high 
concentrations of calcium, conductivity and hardness, indicated near the 
southern corner of the leaching bed, could probably be due to the deposits 
of limestone in that part of the site. 

4.1.4 Concentration Variation with Time 

Figures 21, 22, 23, 24, 25, 26 and 27 indicate contaminants 
concentration variation with time in relation to the natural background 
levels and feed concentrations for the loading period, at the sampling 
well locations where highest concentration of contaminants were recorded. 

4.1.5 Concentration Variation with Distance 

Figures 28, 29, 30, 31, 32, 33 and 34 show the concentration 
variations from C/Cq = 1.0 to 0.1, with respect to distance from the 
centre of the leaching bed, for selected weeks. In general, concentra- 
tion of the contaminants decreased with an increase in distance away from 
the leaching bed. 

The maximum displacement of contaminants in week 145 at the 
concentration level C/C^ = 0.1 was 85 m for chloride, 82 m for nitrate, 
96 m for sulphate, 73 m for sodium, 104 m for hardness, 98 m for calcium 
and conductivity. 

4.1.5 Concentration Variation with Depth 

Figures 35, 36, 37, 38, 39, 40 and 41 show the concentration 
levels of the contaminants at selected depths. Each concentration value 
represents an average of some 40-46 samples collected at each sampling 
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TABLE 16 



MAXIMUM CONCENTRATION LEVELS OF CONTAMINANTS RECORDED IN THE 
SAMPLING WELLS AFTER SUBSURFACE DISPOSAL 



Concentration of Contaminants in mg/L 



Well 
No. 



Nitrate 



Chloride 



Sodium 



Calcium 



Sulphate 



Conduc- 
tivity '' 



Hardness 



t 


0.2 


4 


a 


30 


8 


176 


95 


2 


4.4 


6 


s 


39 


30 


248 


122 


S 


0.7 


7 


i 


88 


22 


490 


243 


4 


5.7 


24 


1 


72 


67 


490 


188 


i 


0,7 


7 


i 


44 


20 


290 


138 


i 


0,5 


6 


4 


41 


26 


261 


134 


7 


0.8 


5 


4 


46 


20 


253 


122 


i 


1.2 


8 


4 


39 


24 


245 


120 


9 


1,6 


8 


7 


86 


30 


500 


257 


10 


8.3 


76 


11 


146 


60 


860 


435 


11 


15.0 


102 


41 


212 


90 


1070 


532 


12 


10.0 


45 


5 


92 


34 


580 


288 


13 


5.4 


31 


1 


91 


33 


490 


296 


14 


2.7 


11 


4 


50 


24 


320 


158 


15 


4.2 


6 


f 


54 


28 


394 


174 


16 


1.4 


19 


7 


115 


41 


650 


345 


17 


6.4 


39 


15 


118 


57 


710 


372 


18 


17.0 


96 


54 


130 


72 


900 


400 


19 


24.0 


100 


69 


141 


70 


970 


414 


20 


21.0 


96 


53 


154 


69 


1000 


430 


21 


3.4 


54 


8 


132 


37 


700 


383 


22 


13.0 


15 


11 


137 


30 


780 


425 


23 


14.0 


91 


23 


150 


65 


939 


481 


24 


16.0 


74 


11 


122 


57 


730 


364 


25 


27.0 


110 


82 


155 


80 


1060 


466 


26 


22.0 


97 


73 


136 


70 


1060 


396 


27 


20.0 


104 


60 


162 


66 


1110 


480 


28 


13.0 


88 


13 


162 


69 


970 


491 


29 


16,0 


85 


38 


149 


63 


960 


460 


30 


18.0 


87 


20 


148 


63 


910 


429 


31 


22.0 


88 


45 


180 


70 


1080 


517 


32 


21.0 


75 


23 


122 


59 


760 


354 


33 


17,0 


75 


6 


100 


54 


710 


348 


34 


1.2 


7 


7 


60 


30 


380 


199 


35 


0.4 


7 


S 


94 


46 


510 


286 


36 


0.6 


6 


i 


104 


28 


533 


296 



Highest 27 110 82 

Value 

* Conductivity reported in ymhos/cm 



212 



90 



1110 



532 
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FIGURE 35 CONCENTRATION OF NITRATE AT SELECTED DEPTHS 

BELOW THE WATER TABLE <mg/L) 
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FIGURE 36 CONCENTRATION OF CHLORIDE AT SELECTED DEPTHS 

BELOW THE WATER TABLE Cing/L> 
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FIGURE 37 CONCENTRATION OF SODIUM AT SELECTED DEPTHS 

BELOW THE WATER TABLE <mg/L> 
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FIGURE 38 CONCENTRATION OF CALCIUM AT SELECTED DEPTHS 

BELOW THE WATER TABLE Cmg/L> 
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FIGURE 39 CONCENTRATION OF SULPHATE AT SELECTED DEPTHS 

BELOW THE WATER TABLE Cmg/L> 
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FIGURE 40 CONCENTRATION OF HARDNESS AT SELECTED DEPTHS 

BELOW THE WATER TABLE <mg/L) 
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FIGURE 41 CONCENTRATION OF CONDUCTIVITY AT SELECTED 
DEPTHS BELOW THE WATER TABLE <p mKo*) 
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level of a multi-level well, between May 1977 and April 1979. The 
figures show two cross-sections x and y taken at the centre of the 
multi-level wells in bays A, B, D and G, F, E, respectively, parallel 
to the length of the leaching bed. 

In general, the highest concentration of contaminants was 
found within one metre of depth below the groundwater table. The level 
of concentration decreased with increased depth below the water table. 
Some of the figures show a variation to the above statement in cases of 
concentrations above the original water table line. This could be due 
to the fluctuations of the water table between the original and maximum 
level of the groundwater, presence of clay lenses and changes in the 
density of the soil medium. In the case of calcium and hardness, the 
increase in concentration with depth could probably be due to ion 
exchange with sodium and potassium on the surfaces of clay particles. 
The decrease in concentration with depth of sodium could be due to ion 
exchange as well as dispersion and dilution. 

4.2 Movement of Contaminants in Groundwater 

The movement of contaminants in groundwater depends on various 
factors, such as, hydraulic loading, subsoil conditions, direction of 
groundwater flow and hydrogeological factors, etc. The travel and 
distribution patterns of the contaminants were the net result of the 
above-mentioned factors in each case. 

'+.2.1 Breakthrough Field Velocities 

Breakthrough at a sample location occurs when the concentration 
level of the element or the dissolved salt in the groundwater increases 
above the background value at that location. The velocity is the rate 
at which a substance travels between two points in unit time. The 
equation used in this study for plotting contours and dispersion plumes 
was : 
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C - C^ 
Concentration = _ = <1 
^o" ^b 

where C is contaminant concentration in the groundwater (after 

disposal) 

C, is background or natural concentration of contaminant 

in the groundwater 

C is contaminant concentration in the feed 
o 

1 is taken as complete saturation 

In this study, breakthrough is taken as the 0.1 concentra- 
tion level. Breakthrough contaminant concentrations were plotted 
at 10 week intervals and velocities were calculated in relation to 
distance travelled in metres per day. 

The breakthrough time and velocity were calculated from the 
contaminant front plotted in accordance with the data collected from 
the sampling well network , with regard to concentration variations 
from C/C = 1.0 to 0.1 of the contaminants, with respect to time and 
distance, as shown in Figures 42, 43, 44, 45, 46, 47, 48 and 49. 

Table 17 shows the breakthrough field velocities of chloride, 
nitrate, sodium, calcium, sulphate, conductivity and hardness for a 
10 week time interval. 

4.2.2 Adsorption Capacity of Soil 

Adsorption may be defined as the tendency exhibited by all 
solids to condense upon their surfaces a layer of any solute or a gas 
with which they are in contact. 

Adsorption capacity of the soil, as used in this study, 
refers to the selective uptake and storage of contaminants (chemicals) 
present in the soil water, per unit weight of the soil. 
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SULPHATE 




FIGURE 46 



FIGURE 47 



CONTAMINANT BREAKTHROUGH FRONT AT 10 WEEK INTERVALS 
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TABLE 17 
BREAKTHROUGH FIELD VELOCITIES OF THE SELECTED CONTAMINANTS 



Contaminants 


Chloride 


Nitrate 


Sodium 


Calcium 


Sulphate 


Conductivity 


Hardness 


Time in 


Weeks 


Velocity 
(m/d) 


Velocity 
(m/d) 


Velocity 
(m/d) 


Velocity 
(m/d) 


Velocity 
(m/d) 


Velocity 
(m/d) 


Velocity 
(m/d) 


10 




0.100 


0.117 


- 


0.104 


,. 


0.113 


0.108 


m 




0.155 


0.152 


- 


0.201 


0.033 


0.171 


0.159 


m 




0.162 


0.157 


O.Oi+2 


0.186 


0.089 


0.179 


0.158 ' 


m 




0.163 


0.160 


0.076 


0.169 


0.095 


0.172 


0.155 


m 




0.158 


0.156 


0.086 


0.155 


0.096 


0.157 


0.149 


m 




O.li+7 


0.11+5 


0.094 


0.141 


0.103 


0.145 


0.138 


m 




0.140 


0.135 


0.098 


0.132 


0.102 


0.134 


0.128 


m 




0.132 


0.128 


0.096 


0.123 


0.100 


0.125 


0.121 


iO 




0.123 


0.122 


0.097 


0.116 


0.099 


0.118 


0.114 


100 




0.115 


0.117 


0.096 


0.111 


0.097 


0.112 


0.109 



Average Velocity 0.140 



0.139 



0.086 



0.144 



0.090 



0.143 



0.134 



- 9U - 



The adsorption capacity of soil could be calculated using an 
extension of the Langmuir equation (4, 14, 35, 38). The modified form 
of the equation is given as follows :- 

V Ap^ 

/ = i^rr (^> 

X o 

or A =(^-l)x — xC^ (b) 

X b 

where 'A' is the adsorption capacity of the soil in 

micrograms per gram 

'V ' is the velocity of the waterfront 

'Vj^' is the velocity of the contaminant front 

'f ' is the porosity of the soil 



'p ' is the bulk density of the soil in grams 



per cubic centimetre 



'C ' is the equilibrium concentration of contaminant 
discharge to the ground, in milligrams per litre 



In addition to the inherent assumptions in the derivation of 
Langmuir 's equation, (4, 14, 35, 38), the following assumptions were 
made for calculation of adsorption capacity of the soil : 

1. The chloride front (case 1) and conductivity front 
(case 2) were not affected by the soil. Hence, the 
adsorption capacity of the soil for chloride (case 1) 
and conductivity (case 2) = 0. 

2. The velocity of the waterfront is equal to the velocity 

of the chloride front (V = V , for case 1) and velocity 

w cl 

of the conductivity front (V = V for case 2), 
■^ w con 

3. The average breakthrough velocity of the contaminant is 
taken as the velocity of the contaminant front (V ). 
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^* The mean concentration value of the contaminant in the 
feed is taken as the equilibrium concentration of the 
contaminant, discharged to the ground (Cq). 
The estimated value of the bulk density (Pj^) for the project site was 
1.6 grams per cubic centimetre and the average bulk porosity (f) as 
measured from the saturated soil samples was 0.35. 

The calculated velocity of the calcium front was highest most 
probably as a result of the desorption of the calcium from limestone 
deposits in the project area. After the calcium, the highest calculated 
value of the contaminant velocity front was that of conductivity. For 
the sake of comparison, two values of adsorption capacity of the soil in 
the project area were calculated based on the assumption that the soil 
adsorption capacity for chloride (case 1) and conductivity (case 2) was 
zero. 

Case 1 : When the adsorption capacity of the soil for chloride 

was assumed equal to zero, the velocity of the waterfront 'V ' was equal 

w ^ 

to the velocity of the chloride front. The value of the calculated 
adsorption capacity 'a' in yg/g was for chloride, 0.03 for nitrate, 
3.21 for hardness, 8.15 for sulphate, and 10.75 for sodium. 

Case 2 : When the adsorption capacity of the soil for conduc- 
tivity was assumed equal to zero, the velocity of the waterfront 'V ' 

w 

was equal to the velocity of the conductivity from 'V ' . The value of 

con 

the calculated adsorption capacity 'A' in pg/g was 0.49 for chloride, 
0.14 for nitrate, 4.81 for hardness, 8.64 for sulphate, and 11.34 for 
sodium. 

The calculated adsorption capacities of the selected contaminants 
based on the above-mentioned two assumptions are presented in Table 18. 
The difference between these two sets of adsorption capacity values is 



TABLE 18 
ADSORPTION CAPACITY OF NORWOOD SOIL FOR THE SELECTED CONTAMINANTS 



Contaminants 


Q 

A vg . Con cen t ra - 
tion in Feed 
(ing/L) 


V ** 

w 

Velocity of 
Front (m 


Water 

/d) 

V 
con 


V *ftft 

X 

Velocity of Con- 
taminant Front 
(m/d) 


A 
Adsorption Capacity 

(yg/g) 






V = V , 
w cl 


\ 


= V 
con 




Chloride (as Cl) 


99.7 


0.140 




0.143 




0.140 









0.49 




Nitrate (as N) 


20.8 


0.140 




0.143 




0.139 




0.03 




0.14 




Sulphate (as SO^) 


63.8 


0.140 




0.143 




0.090 




8.15 




8.64 




Sodium (as Na) 


74.4 


0.140 




0.143 




0.086 




10.75 




11.34 


1 


Calcium (as Ca) 


112.5 


0.140 




0.143 




0.144 




-ve*^*'**!^ 




_ve>'<^"'"'f 


Ol 


Hardness (as CaCo3) 


311.5 


0.140 




0.143 




0.134 




3.21 




4.81 




Conductivity 


974.6 


0.140 




0.143 




0.143 




^ 




Q 





* Mean value is used for Nitrate and Median value used for all other contaminants 
*" Velocity of the waterfront is assumed to be approximately that of chloride and conductivity 
**ft Breakthrough velocities of the contaminants are used 
**** Negative value indicates desorption of containment. 
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negligible. The adsorption by the soil of chlorides, nitrates, 
hardness and conductivity is not very significant. The adsorption of 
sulphates and sodium is higher compared to chlorides, nitrates, hardness 
and conductivity values. In the case of calcium, a desorption is 
indicated. 

'+.2.3 Travel of Contaminants 

Figures 50, 51, 52, 53, 54, 55 and 56 indicate the distance 

travelled by the selected contaminants at the concentration C/C = 0.1, 

o 

with respect to time. Figure 57 shows the relative travel of all the 
selected contaminants. 

Initially, except in the case of sodium, the distance of 
contaminant travel from the edge of the leaching bed increased with time. 
Later, the rate of increase in the distance gradually decreased until 
there was no further increase. The curve shown in the above figures 
became asymptotic and the distance of the contaminant travel at the con- 
centration C/C^ = 0.1 became constant for the rest of the period. 

Figures 50 and 57 initially indicate an increase in distance 

travelled by nitrate with time, at the concentration level C/C = 0.1. 

o 

After reaching a distance of 82 m, the curve became asymptotic. 

Figures 52 and 57 indicate an initial adsorption and then 
desorption of sodium until the curve becomes asymptotic at a distance of 
73 m. 

Figures 53, 55 and 57 indicate an increased travel distance 
of calcium and hardness due to desorption. 

4.2.4 Spread of Contaminants 

The concentration plumes of the selected contaminants were 
plotted, using the concentration levels from C/C^ = 0.9 to 0.1, with 
respect to time. In the case of Figures 62c, 52d, 63a, 63b, 63c, 63d 
and 6Ud, the contours near the edges of the diagram are artificially 
compressed as a result of boundary conditions of the computer plotter. 
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FIGURE 51 TRAVEL OF CHLORIDE AT CONCENTRATION C/Co-0. 1 
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FIGURE 53 TRAVEL OF CALCIUM AT CONCENTRATION C/Co-0. 1 
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FIGURE 54 TRAVEL OF SULPHATE AT CONCENTRATION C/Co-0. 1 
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FIGURE 55 TRAVEL OF HARDNESS AT CONCENTRATION C/Co=»0. 1 
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FIGURE 56 TRAVEL OF CONDUCTIVITY AT CONCENTRATION C/Co-0. 1 
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chloride 

Figures 58a,, 58b, 58c and 58d indicate the distributioii 
pattern of chloride for weeks 17, 63, 97 and 146, respectively,, at 
¥arious ooncentration levels. The chloride plume in the 17th week 
had its peak adjacent tO' well 25 and seemed to move in the southern 
direction. In week 63, the plume expanded in all directions and more 
in the southern direction. In week 97, the spread became more uniform 
and showed movement toward so^uthern and western directions. In week 
146, the spread was more in the northwestern direction, towards the 
river. 

The breakthrough velocity of chloride at the concentration 
C/Cq = 0.1 ranged from 0.10 m/d to 0,16 m/d. The maximum distance 
travelled from the centre of leaching bed was 85 m. The maximumi area 
of chloride spread in week 146 was 2.. 3 hectares. 

Nitrate- 

Figures 59a, 59b, 59c and 59d indicate the movement pattern 
of nitrate in weeks 17, 63, 116 and 146, respectively. The nitrate 
plume in week 17 had its peak of concentration near well 25 and the 
direction of its movement was in the sO'Uthem and southeaste'ro diree-' 
tions. In week 63, the concentration pe^ak was between wells 19 and 20. 
The spread covered a larger area and the direction of move was more in 
the southeastern and northeastern directions. In weeks 116 and 146, 
the spread of concentration contour C/C^ = 0.1 was more in the south- 
western and northwestern directions. 

The breakthrough velocity of nitrate ran,ged from 0.12 m,/d 
to 0.16 ra/d. The maximum distance travelled was 82 m. The maximum 
area of nitrate spread in week 146 was 2.1 hectares. 
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Sulphate 

Figures 60'a„ 60b, 60c and 60d indicate the distribution 
pattern of sulphate in weeks 19 » 683 85 and l^+B „ respectively ,, with 
the peak of the concentration near well 25 in week 19, The concentra- 
tion contour C/C = 0.1 of the sulphate plume seemed tO' mO'Ve in south-- 
eastern and northeasterB directions. In week 63, the movement was 
more in the soutlieastern and northwestern directions. In weeks 85 and 
146, the distribution was more in the southern and western directions. 

The bre.ak through velocity of sulphate ranged from 0'.03 m/d 
to '0.10 m/d. The raaximum distance travelled was 96 m. The maximum 
area of sulphate .spread in well 14-6 was 2.9 hectares. 

Sodium 

Figures 61a 5 61b, 61g and 61d indicate the movement pattern 
o^f sodium, in weeks 53, 85, 116 and 146, respectively. The concentration 
pes-ak was near well 25 in week 53 and the concentration contour C/Cq, = O'.l 
seemed tO' move more in the southe'm direction. The plume expanded in 
all directions in weeks 85 and 116, and showed more spread in the western 
and southern directions in week 146. 

The breakthrough velocity of sodiumi ranged from, 0.0i4 m/d to 
O'.IO m/d. The maximum distance travelled was 73 m,. The maximum area of 
sodium spread in week 146 was 1.7 hectares. 

Calcium 

Figures 62a, 62b, 62c and 62d show the distribution pattern- 
of calcium in weeks 19,, 63, 116 an,d 146. Four concentration peaks 
indicated in week 19, formed one peak, around wells 19, 20', 25 and 26, 
in week 63. The plume seemed to spread in all directions -and moved in 
the southern direction,, The plume eiqpanded in weeks 116 and 146. In 
week 146, more movement was indicated in the southe'm direction. 



- 105 - 



WEEK 19 



WEEK 63 





> ♦ 



FIGURE 60a 



FIGURE 60b 



SULPHATE CONCENTRATION CONTOURS 




SCALE 



WEEK 146 




FIGURE 60O 



METRES 00 



FIGURE 60d 



- 106 - 



WEEK 53 




• * V 



0.1 



* * * ♦ 



• * * * 



WEEK 85 



gill ' 



FtH : 




* • * * 



FIGURE 61a 



FIGURE 61b 



SODIUM CONCENTRATION CONTOURS 



WEEK 116 






h^ I I 




• •B.l • • 



* * # * 



SCALE 



WEEK 146 




FIGURE 61o 



METRES 



-ia 



FIGURE 61d 



- 107 - 



WEEK 19 





FIGURE 62a 



FIGURE 62b 



CALCIUM CONCENTRATION CONTOURS 




WEEK 146 



FIGURE 62o 



SCALE 



METRES 




00 



FIGURE 62d 



- 108 - 

The breakthrough velocity of calcium ranged from 0.10 m/d to 
0.20 m/d. The maximum distance travelled was 98 m. The maximum area 
of calcium spread in week 146 was 3.0 hectares. 

Hardness 

Figures 6 3a, 63b, 63c and 63d indicate the distribution pattern 
of hardness in weeks 19, 70, 98 and lM-6, respectively. Multiple concen- 
tration peaks indicated in week 19 formed one peak and the plume seemed 
to expand more in the southerly direction in weeks 70 and 98. In week 
146, expansion of plume was in all directions and more dispersion in the 
southerly direction was indicated. 

The breakthrough velocity of hardness ranged from 0.11 m/d to 
0.16 m/d. The maximum distance travelled was lOU m. The maximum area 
of spread in week 1U6 was 3.4 hectares. 

Conductivity 

Figures 64a, 64b, 64c and 64d indicate the movement patterns 
of conductivity in weeks 19, 55, 116 and 146, respectively. Two peaks 
of conductivity values in week 19 were located near wells 19 and 25. 
The dispersion was more in the southern direction. In week 55, the peak 
of the concentration was around well 25 and the spread was more in the 
northwestern and southern directions. In weeks 116 and 146, the plume 
expanded in all directions and more in the southern direction. 

The breakthrough velocity of conductivity ranged from 0.11 m/d 
to 0.18 m/d. The maximum distance travelled was 98 m. The maximum area 
of conductivity spread in week 146 was 3.0 hectares. 
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4. 3 Discussion on Groundwater Contaniination 

The groundwater contamination in the study area, due to 
subsurface disposal of sewage effluent from a large size leaching bed, 
may be considered from two main points: magnitude and extent. 

1 . Magnitude of the Problem 

Comparing the concentrations of the selected contaminants 
before and after subsurface disposal, it could be seen that the maximum 
increase in the level of the selected contaminants concentration in the 
groundwater was approximately within one metre below the water table, 
under and in the vicinity of the leaching bed. It may be assumed that 
with time, the concentrations in the vicinity of the bed would eventual- 
ly become equal to the concentrations of the chemicals in the feed. 
Hence, the groundwater in this area is not safe for human consumption 
due to the high level of the contaminants concentration. 

2. Extent of the Problem 

Considering the movement of the selected contaminants at the 
concentration C/Cq = 0.1, the maximum distance travelled and the area 
of spread of the contaminants by week 14-6, were as shown in the Table 19. 

The velocities of the selected contaminants, over the period 
of the study, ranged from O.M- - 0.2 m/d, due to the many factors 
that affected the movement of contaminants in the groundwater, such as, 
the fluctuations in the level of groundwater due to mounding as a 
result of hydraulic loading, precipitation and evapotranspiration losses, 
soil characteristics, presence of rock or clay lenses, attenutation 
reactions, etc. 

Tn accordance with Darcy's law, the velocity of the contami- 
nants increased with the rise in the groundwater level due to mounding. 
With an increase in distance from the peak of the mound, the velocities 
tended to decrease with time. 
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TABLE 19 
SPREAD OF THE SELECTED CONTAMINANTS AT THE 
CONCENTRATION C/Cq =0.1 



Breakthrough Field Maximuni Distance Area of 
Contaminants/ Velocity-'' Travelled from the Spread in 
Parameters (m/d) Centre of the Bed (m) Hectares 



12 2.1 

m^ 2.3 

71 1.7 

ft 3.0 

i§: 2 . 9 

104 3.4 

fS 3.0 



Nitrate 




0.139 


Chloride 




0.140 


Sodium 




0.086 


Calcium 




0.144 


Sulphate 




0.090 


Hardness 




0.134 


Con duct ivi 


■ty 


0.143 



5'fAverage Value 



The average Breakthrough Field velocities of the contaminants 
are indicated in Table 19. For the period of 134 weeks (938 days), 
after subtracting the initial 12 week time lapse before the actual sub- 
surface disposal of sewage effluent began , the maximum radius of the 
contaminants spread from the centre of the leaching bed at the concentra- 
tion level C/Cq =0.1 was calculated to be 104 m. The maximum area of 

2 

the contaminants spread at the concentration C/C^ = 0.1 was 33,975 m or 

3.4 ha as indicated in Table 19. 

Besides the hydraulic loading and hydrogeological factors (9), 
the concentration of the contaminants reaching the groundwater is a 
function of the quality and quantity of the effluent discharge (1, 22, 
34). If the quality of the effluent disposed deteriorates or the 
quantity of the effluent increases in relation to groundwater, the 
contamination would spread over a larger area. 



I"-"! 
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Installation of large subsurface sewage disposal systems , 
hence , should only be allowed in aquifers of such soil and hydrogeologic 
conditions that will provide sufficient attenuation and/or dilution to 
the contaminants, for protection to public health, especially from 
nitrate contamination (2, 20, 26). 



SUMMARY 



- 114 - 

5 . SUMMARY 

An extensive field investigation was conducted to determine 
the environmental impact of a large subsurface sewage disposal system 
on groundwater mounding and contamination. 

The data concerning the fluctuation of the water table beneath 
a large leaching bed were obtained for a period of approximately 2^ 
years from August 1976 to April 1979. The mounding of the water table 
was influenced by the soil conditions, the hydraulic loading rate, 
natural environmental factors, such as, precipitation, evapotranspira- 
tion, and the quantity of infiltration attributed to the spring thaw. 

At the high loading rate of 122,700 L/d, the groundwater 
mound reached 1.99 m in 35 days. For the hydraulic loading rate of 
40,900 L/day, the height of the mound was 1.54 m averaged over 26 weeks 
during the summer and fall seasons in 1978. At the spring thaw period 
in 1978, the maximum rise of the mound as a result of both the hydraulic 
loading and the spring thaw infiltration was 2.75 m. 

Two mathematical models, (Hantush's and Sykes'), were used to 
make predictions of the groundwater mound. Reasonable comparison of 
experimental data and theoretical results were obtained for the follow- 
ing values of input parameters : 

-4 
k - 5 X 10 cm/s, f = 0.2 and H = 12 m. 

From the mathematical modelling exercise, some design guide- 
lines were suggested for the prediction of the groundwater mounding 
beneath a leaching bed. 

In order to study the groundwater contamination due to sub- 
surface sewage effluent disposal from a large leaching bed, seven 
chemical parameters, namely: nitrate, chloride, sodium, calcium, 
sulphate, hardness and conductivity were selected. 
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The concentration of contaminants under and in the vicinity 
of the leaching bed increased significantly after effluent disposal. 

The highest contaminant concentration with respect to depth 
was found within one metre below the groundwater table. A decrease in 
concentration levels with depth was indicated. 

Contaminants concentration tended to decrease with an increase 

in distance away from the leaching bed. The maximum distance travelled 

at week 146, at the concentration C/C = 0.1 was 85 m for chlorides, 

o ' 

82 m for nitrates, 96 m for sulphates, 73 m for sodium, 98 m for calcium, 
104 m for hardness and 98 m for conductivity. The distance and time 
curves of the contaminants became asymptotic after reaching the above 
distances. 

The velocity of the contaminants at the concentration C/C = 0.1 
varied from O.IU m/d for chloride and nitrate, 0.09 m/d for sulphate and 
sodium, 0.14 m/d for calcium and conductivity, and 0.13 m/d for hardness. 

A modified version of the Langmuir Equation was used to deter- 
mine the adsorption capacity of the soil. The adsorption capacity of the 
soil in yg/g, for the case V = V , was for chlorides, 0.03 for 
nitrates, 3.21 for hardness, 8.15 for sulphates and 10.75 for sodium. The 
adsorption capacity of the soil in yg/g, for the case V^ = V , was 0.49 
for chlorides, 0.14 for nitrates, 4.81 for hardness, 8.64 for sulphates 
and 11.34 for sodium. 

In general, the contaminants seemed to spread in all directions 
in the beginning and more distribution was indicated toward the river in 
the later stages. 
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6.0 CONCLUSIONS 

A. Conclusions related to the groundwater mounding portion of 
the study are as follows :- 

1. Many factors affected the fluctuations of the water table 
beneath a large leaching bed including precipitation, evapotranspira- 
tion, and the hydraulic loading rate. The relative significance of 
these factors on the movement of the water table depends on the site 
and on the seasons. In the case of Norwood, the influence of the 
infiltration at the spring thaw period in 1978 was more significant 
than the effect of the hydraulic loading (40,900 L/d) . 

2. Two mathematical models (Hantush's and Sykes'), were used to 
predict the extent of the groundwater mounding. For identical values 
of input parameters similar results were given by the two models. 
Because of the relative simplicity in the use of the Hantush method, it 
is more suitable as a routine design tool. However, the flexibility of 
Sykes' finite -element model should not be overlooked, and the numerical 
method may be used for more complicated site conditions if the soil 
characteristics can be measured with greater accuracy. 

i« In using Hantush's solution, it is recommended that the 
'falling head' and the 'recovery' tests be used to measure the 'k' 
(permeability) values of sandy silt soils. In addition to the measure- 
ments of 'k' and specific yield (s), at individual locations, a detailed 
study of the soil stratigraphy at the proposed site should be carried 
out so that a proper interpretation of the soil test data can be made. 

%. The Hantush method can be used to predict the mounding due 
to a constant hydraulic loading and it cannot be used to predict the 
mounding caused by the spring thaw infiltration. From this study, it 
is concluded that an estimate of the rise in water table due to the 
infiltration of surface water can be made on a statistical basis 
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according to the weather data and a ininiinuiii of one year's observational 
data on the natural fluctuations of the water table. The magnitude of 
the maximum water table rise should be equal to the sum of the rise 
caused by the hydraulic loading and due to surface infiltration. 

B. Conclusions related to the groundwater contamination section 
of the study are as follows : - 

1 . Comparing the concentration of chemicals in the groundwater 
before and after subsurface disposal, the groundwater under and in the 
vicinity of leaching bed approached the chemical concentrations of the 
feed. 

2. The maximum concentration of the contaminants with respect to 
depth was within one metre below the water table. The concentrations tend 
to decrease with an increase in depth. 

3. Groundwater mounding increased the contaminant velocity due 
to increased hydraulic gradient. 

4. The maximum distance travelled by the contaminants at the 
concentration C/C =0.1 observed for 13^ weeks was 104 m from the 
centre of the leaching bed; the maximum area of the contaminants spread 
was approximately 3.4- ha. 

5. The time and distance curves plotted for the selected contami- 
nants at the concentration level C/C = 0.1 became asymptotic after 
reaching the maximum distances as indicated in week 14-6. 

6. The magnitude and extent of the groundwater contamination 
depended on several factors such as the hydraulic loading rate, quantity 
and quality of sewage disposed, climatic and hydrogeological factors, 
subsoil conditions and size of the aquifer. 
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7.0 RECOMMENDATIONS 

Based on the above conclusions, it is recommended that: 

1. Large subsurface sewage disposal systems should only be 
installed in aquifers of such soil and hydrogeologic conditions that 
will provide sufficient attenuation and/or dilution to the ground- 
water contaminants to keep their concentration levels within safe 
federal and provincial limits. 

i# These systems should be designed with due regard to the 
hydrogeological factors of the site (as mentioned in this report) 
so that hydraulic failure due to groundwater mounding does not occur, 
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1 . INTRODUCTION 

To model the secondary effluent discharge from a leaching 
bed located in the unsaturated zone and the groundwater mounding that 
results, the governing flow equation for a three dimensional spatial 
domain is: 

(C + i s ) 1^ = ^ [k. . (0) II- + k._ (6)] i, j = 1, 3 (1) 

n s dt dx. ij ax. 13 

where C is the specific moisture capacity, a function of the moisture 
content 9 and the pressure head i{j; n is the porosity; S is the 
specific storage; and k.. (0) is the hydraulic conductivity tensor given 
here as a function of the moisture content. For the saturated zone, 
k. . (0) yields the saturated hydraulic conductivities. Equation (1) is 
nonlinear, in that the parameters C and k.. (9) are functions of the 
solution set. The specific storage coefficient Sg accounts for the 
compressibility of both the pore fluid and the formation itself with 
vertical consolidation being assumed for the latter. 

The use of equation (1) to model the flow system results in 
the discharge from the leaching bed being modelled via either the 
boundary conditions associated with the equation or point source terms. 
Infiltration and evapotranspiration at the air-soil interface is also 
included in the model through boundary conditions. 

To solve equation (1) and the associated initial conditions 
and boundary conditions, numerical solution techniques are required 
with the finite difference method and the finite element methods being 
commonly used. (Freeze, 1971; Dillon et. al. , 1978; Segol, 1976). 
For three dimensional analyses, such models are both difficult and 
expensive to use. In addition, it is often the case that parameter 
uncertainty and a lack of complete spatial definition results in the 
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fact that three dimensional phenomena can often be modelled in two 
dimensions with little or no loss of information. For systems in 
which vertical effects do not predominate, areal models can be 
developed for the saturated zone by suitably integrating equation (1) 
over the vertical domain yielding an equation of the form: 

where T^. = K^."b is the transmissivity tensor; K. . is the saturated 
hydraulic conductivity tensor averaged over the aquifer thickness b; 
h is the head; S is the storativity and Q(x. ) is a source term. The 
boundary conditions associated with the equation are: 

h(x^) = h(x^) on r^ (3a) 

and q(x^) = q(x^) on T^ (3b) 

where h(x. ) is a prescribed head on boundary V and q(x. ) is a 
prescribed flux on boundary T^ where T + ^ ^ = T is the external 
boundary to the modelled domain. In the use of equation (2), the 
source/sink term Q(x. ) can vary spatially and can represent infiltra- 
tion, evapotranspiration or recharge loading from a leaching bed. For 
mounding of phreatic aquifers, the aquifer thickness as modelled by b 
is a measure of the height between the water table and the underlying 
confining bed. As such, it is a function of the solution set h(x. , t). 

Equation (2) and its governing boundary and initial conditions 
can be solved both exactly and numerically. For the former, simplify- 
ing assumptions are commonly made and, in fact, are often necessary. 
Generally, the transmissivity is assumed to be constant. In addition, 
solutions are often available for regular spatial domains only. 
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2* SYKES FINITE ELEMENT MODEL 

As is the case with the three-dimensional equation (1), the 
numerical solution of equation (2) is generally based on either the 
finite difference or the finite element methods. The literature 
abounds with such models and in this study, a finite element model 
developed by Sykes (1975) is used. The model is highly flexible and 
can investigate contaminant transport as well as flow in either two- 
dimensional variably saturated cross-sections. A wide range of 
boundary conditions is available in the model. 

To solve equation (2), the Galerkin method of weighted 
residuals is used to replace the given differential equation by a set 
of integral equations which are solved numerically. In this sub- 
section, the Galerkin equation for fluid flow in two-dimensional 
areal spatial domains is developed. Details of the finite element 
method used to solve the equation are also presented. 
Galerkin Flow Equation 

In the Galerkin procedure , the general approach is to assume 
a trial solution whose functional dependence on position is chosen, 
but which includes undetermined functions of time. For equation (2), 
the trial solution assumed has the form: 

h = h® = <n> {h} (4) 

where <n> is a set of prescribed non-dimensional spatial interpolation 
functions satisfying the boundary condition associated with equation 
(2) and h^ is the set of undetermined time dependent heads defined at 
the nodel points at a discretized continuum defining the solution 
domain. The sets <n> and {h} must have the same number of ordered 
members . 
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Differentiation of equation (4) with respect to the spatial co-ordinates 
yields the pressure head gradients. 

9h ^9n . r, -1 , ^ 

1 1 

in which ^"^^^ is the transform vector giving the gradients in terms 

i 
of the set of nodal point variables. 

To develop the Galerkin equation, the flow equation (2) is 

expressed in terms of the partial differential operation L(h) giving 

^ t = al: "i: 1^' - <5(x,) = L(h) (6) 

where i, j = 1, 2. In this equation, S and Q remain constant within an 
element of the subdivided continuum while T.. is expressed as a 
functional coefficient utilizing the same set of non-dimensional inter- 
polation functions as that used in equation (4). 

To determine a solution set at the nodal points at time 
t + At given the solution set at the nodal points of the discretized 
system at time t, a finite difference approximation for the time deriva- 
tive is introduced. Thus, equation (6) is written as: 

h^ '^^ -h^ =aL(h) ^ ^^^ t (1 -a) L(h)^ (7) 



S At 

where the recurrence parameter a is such that o< a ^ 1 . 

The use of the trial solution (4) in equation (6) results in 
a residual 

R(h^) i L(h^) - ^ (8) 
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As the number of members in the set {h} increases, the 
residual will become smaller, with the exact solution being obtained 
when the residual is identically zero. As an approximation to this 
ideal the integrals of the weighted residuals are set equal to zero 
giving : 

/^ {n} R(<n> {h}) dA = (9) 

where following the Galerkin procedure the weights are selected as the 
interpolation function {n}. 

In the manner of equation (9), the Galerkin flow equation 
can thus be developed as: 



At 



ta[B]] {h}^^ ^= h 



[A] + a[B] I {h}" - "^ = I ±_ [A] + (1 -a) [B] 



{h}^ + {c} + {D} 
(10) 



where 



m 
[A] = I / {n} S<n> dA (11a) 

e = 1 



[B] = E / {— } <n> {T } On/9x> 
e = 1 dx 

+ {|E} <n> {T } On/3y> dA (lib) 

3y yy 

m 
{c} = Z / {n} Q (x^) dA (lie) 

e = 1 

(d) = Z ^^ {n} T^^ <||> dy - Cn) T^^ f dx 



^2 



(lid) 



The vector (d) of equation (10) results from the application 
of Green's theorem in a plane to the second order term of equation (6), 
thus reducing it to a first order term and a surface integral. In so 
doing, the Neumann boundary condition of (3b) is incorporated in the 
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Galerkin equation. The nonlinear transmissivities, a function of the 
aquifer thickness b, are evaluated at time t + At/2 and are assumed 
to remain constant over a time step. 
Finite Element Algorithm 

To discretize the spatial domain and to develop the non- 
dimensional interpolation function introduced in equation (4), this 
study uses a mixed isoparametric quadrilateral element. The element 
contains between 4 and 12 degrees of freedom with the sides being 
either linear, quadratic or cubic. The integration of equation (11) 
is accomplished using a Gausslegendre Quadrature. With such a proce- 
dure, the functional coefficients (in this study Ti . ) are evaluated 
at the Gauss points. Two, three or four Gauss points are used in a 
given direction depending on whether the higher order polynomial of 
the two sides in that direction is linear, quadratic or cubic respec- 
tively. 

In the Galerkin flow equation, the matrix { B} is a function 
of the dependent variable h. To obtain an initial estimate of the 
head at time t + At, thus allowing the evaluation of {B}, a linear 
interpolation scheme is used such that for the nodal variable h 

h^ ^ ^^ = h^ ^ (h^ - h^ - ^^) ^ (12) 

At^ 

The Newton Raphson technique (Sykes, 1975) is used to improve 
these estimates of head. In this method, for a nonlinear finite 
element equation of the form: 

[A (h)] {h} - F = (13a) 
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Improvements to the es.tiiiate O'f h " are obtained from'^ 



where 



{h } = {h _,} + [ACh _)r^ (-R) 
new old old 



R = [A(h^,^^)]{h^^^^} = F 



(13b) 
(13c) 



Rapid convergence is achieved by incorporating an Aitken 

accelerator with the Newton-'Raphson technique (Sykes, 1975).. 

-1 

T.O evaluate the matrix form [A(h ^ ,)] (-R) of the 

old 

Newton-Raphson technique ^ an upper banded Gaussian elimination solver 
is used. 

The finite element algorithm is programmed in fortran and 
run on an IBM, 3i 70'. 14-5 machine... 



Information Required for Analysis 
mesh (nodes and elements) 
porosity 

elevation O'f con,fining layer 
permeability (hydraulic conductivity) 
hydraulic loading 
boundary conditions 



